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Abstract. We consider Calderon's inverse problem with partial 
data in dimensions n > 3. If the inaccessible part of the boundary 
satisfies a (conformal) flatness condition in one direction, we show 
that this problem reduces to the invertibility of a broken geodesic 
ray transform. In Euclidean space, sets satisfying the flatness con- 
dition include parts of cylindrical sets, conical sets, and surfaces 
of revolution. We prove local uniqueness in the Calderon prob- 
lem with partial data in admissible geometries, and global unique- 
ness under an additional concavity assumption. This work unifies 
two earlier approaches to this problem f [KSU07] and |Is07j ) and 
extends both. The proofs are based on improved Carleman esti- 
mates with boundary terms, complex geometrical optics solutions 
involving reflected Gaussian beam quasimodes, and invertibility of 
(broken) geodesic ray transforms. This last topic raises questions 
of independent interest in integral geometry. 
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1. Introduction 

This article is concerned with inverse problems where measurements 
are made only on part of the boundary. A typical example is the inverse 
problem of Calderon, where the objective is to determine the electrical 
conductivity of a medium from voltage and current measurements on its 
boundary. The mathematical formulation of this problem is as follows. 
Let Q C M. n , n > 2, be a bounded domain with smooth boundary. 
Given a positive function 7 G L°°(Q) (the electrical conductivity of 
the medium) and two open subsets T^, IV of dfl, consider the partial 
Cauchy data set 

C* D *» = {H rD , 7 cU| r J ; div( 7 V«) = in Q, u G H\Q), 

suppOlan) c r D }. 

The Calderon problem with partial data asks to determine the con- 
ductivity 7 from the knowledge of C^ d,Fn for possibly very small sets 
T^Ttv. Here d u is the normal derivative, and the conormal derivative 
jdviilgn is interpreted in the weak sense as an element of H~ l / 2 (dQ). 

A closely related problem is to determine a potential q G L°°(Q) from 
partial boundary measurements for the Schrodinger equation, given by 
the partial Cauchy data set 

C[°' r " = {(u\ ro , d v u\r N ) ; (-A + q) = in SI, u G H A (Q), 

supp(w| an ) c r D }. 

Here we use the space 

F A (0) = {116 L 2 (n) ; Am G L 2 (n)}, 

and the trace u\qq and normal derivative d v u\oQ are in H^ 1 ^ 2 {dVL) and 
H^ i ^ 2 {dVt) (see [BU01J). Above, one thinks of u\qq as Dirichlet data 
prescribed only on Tjj, and one measures the Neumann data of the 
corresponding solution on T^. If A 7 : H l l 2 {VL) — > H~ l l 2 (dVt) is the 
Dirichlet-to-Neumann map (DN map) given by 

A 7 : u\gn 1 — ^ ■fdvulaa where u G H 1 ^) solves div(7Vw) = in ft,, 

then the partial Cauchy data set is a restriction of the graph of A 7 , 

C r , D ' FN = {(/ir ,A 7 /| r J; / G H^ 2 ^), supp(/) c T D }. 

A similar interpretation is valid for Cg D,rN provided that is not a 
Dirichlet eigenvalue of — A + q in Q. 

The problems above are well studied questions in the theory of in- 
verse problems. The case of full data (T D = T N = dfl) has received 



CALDERON PROBLEM WITH PARTIAL DATA 



3 



the most attention. Major results include |SU87j . |HTllj in dimen- 
sions n > 3 and [Na96j . [AP06] . |Bu08j in the case n = 2. In particu- 
lar, it is known that the set C® n,9n determines uniquely a conductivity 
7 G C 1 (fi) if 72 > 3 and a conductivity 7 G L°°(Q) if n = 2. These re- 
sults are based on the method of complex geometrical optics solutions 
developed in |SU87] for n > 3 and in |Na96j . |Bu08] in the case n = 2. 

The partial data question where the sets Yd or may not be the 
whole boundary has also attracted considerable attention. We mention 
here four approaches, each of which gives a slightly different partial data 
result. Formulated in terms of the Schrodinger problem, it is known 
that C^ d,Vn determines q in Q in the following cases: 

(1) n > 3, the set is possibly very small, and is slightly 
larger than dQ \ T D (Kenig, Sjostrand, and Uhlmann [KSU07] ) 

(2) n > 3 and Y D = T N = T, and dfl \ T is either part of a hyper- 
plane or part of a sphere (Isakov |Is07j ) 

(3) n = 2 and Yd = = T, where T can be an arbitrary open 
subset of dQ (Imanuvilov, Uhlmann, and Yamamoto |IUY10j ) 

(4) n > 2, linearized partial data problem, = Tn = T where 
T can be an arbitrary open subset of dQ (Dos Santos, Kenig, 
Sjostrand, and Uhlmann |DKSjU09|) 

Approaches (l)-(3) also give a partial data result of determining 
7 from C^ d,Fn with the same assumptions on the dimension and the 
sets T D ,T N . In (4), the linearized partial data problem asks to show 
injectivity of the Frechet derivative of A q at q = instead of injectivity 
of the full map q (-)■ A g , when restricted to the sets and Tn- 

It is interesting that although each of the four approaches is based 
on a version of complex geometrical optics solutions, the approaches 
are distinct in the sense that none of the above results is contained in 
any of the others. The result in [KSU07] uses Carleman estimates with 
boundary terms, given for special limiting weights, that allow to control 
the solutions on parts of the boundary, whereas |Is07] is based on the 
full data arguments of |SU87j and a reflection argument. The work 
|IUY10] gives a strong result that only requires Dirichlet and Neumann 
data on any small set, but the method involves complex analysis and 
Carleman weights with critical points and does not obviously extend to 
higher dimensions. Finally, |DKSjU09| is based on analytic microlocal 
analysis but is so far restricted to the linearized problem. 
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Nevertheless, given that there exist several approaches to the same 
problem, one expects that a combination of ideas from different ap- 
proaches might lead to improved partial data results. In this paper 
we unify the Carleman estimate approach of [KSU07J and the reflec- 
tion approach of |Is07j . and in fact we obtain the main results of both 
[KSU07] and [E07] special case. 

The method also allows to improve both approaches. Concerning 
[Is07j, we are able to relax the hypothesis on the inaccessible part 
Tj = dQ \ T of the boundary: instead of requiring Tj to be completely 
flat (or spherical), we can deal with Tj that satisfy a flatness condition 
only in one direction. Compared with [KSU07j we remove the need 
of measurements on certain parts of the boundary that are flat in one 
direction; and in certain cases where dQ may not have any symmetries, 
we eliminate the overlap of Yd and needed in |KSU07] . The method 
eventually boils down to inverting geodesic ray transforms (possibly for 
broken geodesies). In some cases the invertibility of the ray transform 
is known, but in other cases it is not and in these cases we obtain 
a reduction from the Calderon problem with partial data to integral 
geometry problems of independent interest. 

We list here some further references for partial data results, first 
for the case n > 3. The Carleman estimate approach was initiated in 
[BUOlj and [KSU07J. Based on this approach, there are low regularity 
results |Kn06j . |Zhl2j . results for other scalar equations |DKSjU07] , 
[KS07j . [Chl2] and system s [STlO] . stability results |HW06j . and re- 
construction results [NSlOj . The reflection approach was introduced in 
[Is07j, and has been employed for the Maxwell system [COS09J. Partial 
data results for slab geometries are given in [LiUlOj . [KLU12j . Also, 
just before this preprint was submitted, Imanuvilov and Yamamoto 
announced a partial data result for domains Q = Q x (0, 1) stating 
that the Cauchy data set of a potential q £ C a (Q) with Dirichlet and 
Neumann data restricted to T x (0, 1), where T C dQ is an open subset, 
determines the potential in the set (Q \ Qq) x (0, 1) where Qq is the 
convex hull of dfl \ T. This is similar to the results in Section 13.11 in 
this paper. 

In two dimensions, the main partial data result is |IUY10j which 
has been extended to more general equations [IUYllaj . combinations 
of measurements on disjoint sets [IUYllbj . less regular coefficients 
[IY12a] , and some systems |IY12bj . An earlier result is in [ALP05J. 
In the case of Riemann surfaces with boundary, corresponding partial 
data results are given in jGTllaj . jGTllbj . jACTUllj . 

In the case when the conductivity is known near the boundary, the 
partial data problem can be reduced to the full data problem [AU04j . 
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[AK12j . [HPS12J. Also, we remark that in the corresponding prob- 
lem for the wave equation, it has been known for a long time (see 
[KKL01] ) that measuring the Dirichlet and Neumann data of waves on 
an arbitrary open subset of the boundary is sufficient to determine the 
coefficients uniquely up to natural gauge transforms. Recent partial 
results for the case where Dirichlet and Neumann data are measured 
on disjoint sets are in |LO10j . |LQ12j . 

The structure of this paper is as follows. Section [T] is the intro- 
duction. Section [5] states the main partial data results in this paper 
in the setting of Riemannian manifolds, and Section [3] considers some 
consequences for the Calderon problem with partial data in Euclidean 
space. Section H] gives a Carleman estimate that is used to control so- 
lutions on parts of the boundary, and Section [5] discusses a reflection 
approach that can be used as an alternative to Carleman estimates in 
some cases. In Section [6] we give the proofs of the local uniqueness 
results for simple transversal manifolds, based on complex geometrical 
optics solutions involving WKB type quasimodes. In Section [7] we dis- 
cuss a more sophisticated quasimode construction based on reflected 
Gaussian beams, and in Section [5J we show how complex geometrical 
optics solutions involving reflected Gaussian beam quasimodes can be 
used to recover the broken ray transform of a potential from partial 
Cauchy data. 

Acknowledgements. C.K. is partly supported by NSF, and M.S. is 
supported in part by the Academy of Finland and an ERC Starting 
Grant. M.S. would like to thank David Dos Santos Ferreira, Yaroslav 
Kurylev, and Matti Lassas for several helpful discussions and for al- 
lowing us to use arguments from the paper |DKLS12] . which was in 
preparation simultaneously with this manuscript, involving Gaussian 
beam quasimodes and a reduction from the attenuated ray transform 
to the usual ray transform. M.S. expresses his gratitude to the Depart- 
ment of Mathematics of the University of Chicago, where part of this 
work was carried out. 

2. Statement of results 

Our method is based on ideas developed for the anisotropic Calderon 
problem in |DKSaU09] , and even though much of the motivation comes 
from the Calderon problem with partial data in Euclidean domains, it 
is convenient to formulate our main results in the setting of manifolds. 
The Riemannian geometry notation used in this paper is mostly the 
same as in [DK SaU09j . 
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Definition. Let (M, g) be a compact oriented Riemannian manifold 
with C°° boundary, and let n = dim(M) > 3. 

1. We say that (M,g) is conformally transversally anisotropic (or 
CTA) if 

(M,g) cc (R x Mo,g), g = c(e®g ) 

where (Mo,<7o) is some compact (n — l)-dimensional manifold 
with boundary, e is the Euclidean metric on the real line, and 
c is a smooth positive function in the cylinder R x Mq. 

2. We say that (M,g) is admissible if it is CTA and additionally 
the transversal manifold (Mo, go) is simple, meaning that the 
boundary dMo is strictly convex (the second fundamental form 
is positive definite) and for each p e M Q , the exponential map 
exp p is a diffeomorphism from its maximal domain of definition 

in T p Mq onto Mq. 

The uniqueness results in [DKSaU09] were given for admissible man- 
ifolds. In this paper we will give results both for admissible and CTA 
manifolds. In the main results, we will also assume that there is a 
compact (n — l)-dimensional manifold (Mo, go) with smooth boundary 
such that 

(2.1) (M,g) C (R x M ,g) CC (R x M ,g), g = c{e®g ) 

and the following intersection is nonempty: 

DM fl(Ix dMo) ^ 0. 

Under some conditions, it will be possible to ignore boundary measure- 
ments in the set <9Mfl(R x dMo). In the results below we will implicitly 
assume that the various manifolds satisfy ( 12. ip . and if (M,g) is admis- 
sible it is also assumed that (M ,g ) is simple (but (M ,g ) need not 
be simple, since its boundary may not be strictly convex). 

Write x = (xi,x') for points in R x M where X\ is the Euclidean 
coordinate. The approaches of [KSU07] , |DKSaU09] are based on com- 
plex geometrical optics solutions of the form u = e Tip (m + r) where tp 
is a special limiting Carleman weight. We refer to [DKSaU09] for the 
definition and properties of limiting Carleman weights on manifolds. 
For present purposes we only mention that the functions tp(x) = ±Xi 
are natural limiting Carleman weights in the cylinder (R x Mo,g). 

The weight <fi(x) = x\ allows to decompose the boundary dM as the 
disjoint union 

dM = dM + U <9M_ U <9M tan 
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where 

dM± = {xedM; ±d v ip{x) > 0}, 
<9M tan = {x G <9M ; d v tp(x) = 0}. 

Here the normal derivative is understood with respect to the metric 
g. Note that d v tp = on R x dM whenever (M ,g ) CC (M ,g ). 
We think of <9M tan as being flat in one direction (the direction of the 
gradient of tp). For the sake of definiteness, the sets dM± = dM±((p) 
will refer to the weight (p(x) = x\ in this section, but all results remain 
true when dM + and <9M_ are interchanged (this amounts to replacing 
the weight X\ by —X\). 

Next we give the local results for the Calderon problem with partial 
data on manifolds. In these results we say that a unit speed geodesic 
7 : [0, L] — > Mq is nontangential if its endpoints are on 8Mq, the vectors 
7(0), 7(L) are nontangential, and 7(2) G M™ 1 for < t < L. We also 
define the partial Cauchy data set as 

C T g D / N = {Hr D , Mr*) I (~ A 9 + g) = in M, M G H Ag (M), 

supp(M| 9A f) c r D } 

where H Ag (M) = {u e L 2 (M) ; A g u G L 2 (M)} &ndu\ dM G H~ x / 2 (dM), 
d u u\dM ^ H^/ 2 {dM) by the same arguments as in |BU01I . 

To explain the results, it is convenient to think in terms of the fol- 
lowing special case. 

Example. Let M be a compact manifold with boundary consisting of 
three parts, 

M = Mieft U M mid U M right 

where M mic j = [a, b] x M for some compact manifold (Mo, go) with 
boundary, M le f t C {xi < a] x M , and M right C {xi > 6} x M . We 
also assume that 

dM_ = Mi eft n 9M, <9M+ = M right n 9M. 

In this case <9M tan = [a, b] x <9M . 

The methods developed in this paper suggest that it should suffice to 
measure Neumann data on dM + for Dirichlet data supported in dM_, 
with no measurements required on <9M tan . However, in the results 
below we need a part T a C <9M tan that is accessible to measurements, 
and Ti = <9M tan \ T a is the inaccessible part. Suppose for simplicity 
that 

T a = [a,b]xE, T t = [a,b]x(dM \E) 
for some nonempty open subset E of dM . 
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In this setting, Theorem 12.11 implies that from Neumann data mea- 
sured near dM + UT a with Dirichlet data input near dM_ UT a , one can 
determine certain integrals of the potential q in the set 



where the union is over all nontangential geodesies in M with end- 
points on E. Moreover, if the local ray transform is injective in this set 
in a suitable sense, then one can determine the potential in this set by 
Theorem 12.21 Theorem 12.41 shows that one can go beyond this set and 
extract information about integrals of q over all nontangential broken 
rays with endpoints on E, and Theorem 12.31 gives a global uniqueness 
result in the case where dM tan has zero measure. 

Theorem 2.1. Let (M,g) be an admissible manifold as in (12. ip . and 
let gi, q 2 G C(M). Let Tj be a closed subset of <9M tan; and suppose that 
for some nonempty open subset E of OMq one has 



Given any nontangential geodesic 7 : [0, L] — > M with endpoints on 
E, and given any real number \, one has 



Jo 

Here q\ — q 2 is extended by zero outside M, and ( • )" denotes the Fourier 
transform in the x\ variable. 

The previous theorem allows to conclude uniqueness of potentials in 
sets where the local ray transform is injective in the following sense. 

Definition. Let (M ,g ) be a compact oriented manifold with smooth 
boundary, and let O be an open subset of M . We say that the local 
ray transform is injective on O, if any function / 6 C(Mq) with 



RxU#,i]) 



7 





/ dt = for all nontangential geodesies 7 contained in O 



7 



must satisfy f\o = 0. 
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Theorem 2.2. Assume the conditions in Theorem 1 2. 1\ Then qi = q 2 
in M n (M x O) for any open subset O of M such that the local ray 
transform is injective on O and O R dM C E. 

The local ray transform is known to be injective in the next three 
cases (the second case will be used in Section [3]): 

1. (Mo, go) — (^o,e) where Qq C M rt_1 is a bounded domain with 
C°° boundary, e is the Euclidean metric, E is an open subset 
of 8Qq, and O is the intersection of f2 with the union of all 
hyperplanes in IR n_1 that have 8Qq \ E on one side. The com- 
plement of this union is the intersection of half-spaces and thus 
convex. If the integral of / G C(Qq), extended by zero to IR n_1 , 
vanishes over all line segments in O, then the integral over all 
hyperplanes that do not meet dQo \ E also vanishes, and it fol- 
lows from the Helgason support theorem [He99j that the local 
ray transform is injective on O. 

2. (M ,(?o) CC (M ,go) are simple manifolds with real-analytic 
metric, and T is an open set of nontangential geodesies in 
(M ,(7o) such that any curve in T can be deformed to a point 
on 8Mq through curves in T . In such by a result of 
[Kr09] the local ray transform is injective on the set O of all 
points in Mq that lie on some geodesic in J 7 . 

3. If dim(Mo) > 3 and if dM$ is strictly convex at a point p G dMo, 
then p has a neighborhood O in M Q on which the local ray 
transform is injective. This is a very recent result from |UV12j . 

In Theorem I2.2[ if the nontangential geodesies with endpoints on E 
cover a dense subset O of Mq and if the local ray transform is injective 
in O, we obtain a global uniqueness result stating that qi = q 2 in M. 
An example of such a result under a concavity assumption is given in 
Section [3761 The method for proving Theorems 12.11 and 12.21 also gives a 
global result when the set dM tan has zero measure. Then no geometric 
conditions are required and one can recover the potential in all of M. 

Theorem 2.3. Let (M,g) be an admissible manifold and assume that 
qi,q2 G C(M). If dM tan has zero measure in DM, then 



Next we wish to gather information on the potentials beyond the set 
that can be reached by transversal geodesies with endpoints on E. To 
do this, we will use broken geodesies in the transversal manifold that 
go inside M , reflect finitely many times and eventually return to E. 

Definition. Let (M , <?o) be a compact manifold with boundary. 




10 



CARLOS KENIG AND MIKKO SALO 



(a) We call a continuous curve 7 : [a, b] — > M a broken ray if 7 
is obtained by following unit speed geodesies that are reflected 
according to geometrical optics (angle of incidence equals angle 
of reflection) whenever they hit a point of dM . 

(b) A broken ray 7 : [0, L] — > M is called nontangential if 7(t) is 
nontangential whenever ^(t) e OMq, and additionally all points 
of reflection are distinct. 

The next theorem is a generalization of Theorem 12.11 in the sense 
that it allows arbitrary transversal manifolds and recovers integrals over 
all nontangential broken rays (instead of just nontangential geodesies) 
with endpoints on E. However, it is stated with a weaker partial data 
condition. 

Theorem 2.4. Let (M,g) be a CTA manifold as in (12. ip . and let 
q\,qi £ C(M). Let Tj be a closed subset of <9M tan; and suppose that for 
some nonempty open subset E of dM one has 

r.cMx (dM \ E). 

Let Y a = <9M tan \ Ti, and assume that 



where 



rir D ,r N _ f<r D ,r N 
9,91 9,92 



D — J- N 



for some neighborhood T of the set dM + U dM_ U T a in dM. 

Given any nontangential broken ray 7 : [0, L] — > M with endpoints 
on E, and given any real number \, one has 

e- 2Xt (c(q 1 -q 2 )Y(2X, 1 (t))dt = 0. 







Here q\ — qi is extended by zero outside M , and ( ■ )" denotes the Fourier 
transform in the x\ variable. 

It is natural to ask whether a function in Mq is determined by its 
integrals over broken rays with endpoints in some subset E of dM 
(that is, whether the broken ray transform is injective). Combined 
with Theorem 12.41 and with the proof of Theorem 12.2} such a result 
would imply unique recovery of the potential in the whole manifold 
M. However, it seems that there are very few results in this direction, 
except for the case where E is the whole boundary and the question 
reduces to the injectivity of the usual ray transform (see |Sh94] ) . 

Eskin |Es04] has proved injectivity in the case of Euclidean broken 
rays reflecting off several convex obstacles, with E being the boundary 
of a smooth domain enclosing all the obstacles, if the obstacles satisfy 
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additional restrictions (in particular the obstacles must have corner 
points and they cannot be smooth). Ilmavirta |I112j has recently given 
partial injectivity results for the broken ray transform in the Euclidean 
unit disc. See also |FMS11] . |Lol2] for related results. However, the 
following question seems to be open even in convex Euclidean domains 
except when E = OMq. 

Question. Let (M ,g ) be a simple manifold, let E be a nonempty 
open subset of dM , and assume that f e C(M Q ) satisfies 

[ L f( 7 (t))dt = 
Jo 

for all nontangential broken rays 7 : [0, L\ — > Mq with endpoints on E. 
Does this imply that f = ? 

3. The Euclidean case 

In this section, we indicate some consequences of the previous results 
to the Calderon problem with partial data in Euclidean space. We 
assume that Q C 1R 3 is a bounded domain with smooth boundary 
equipped with the Euclidean metric g = e, and q±, q 2 G C(fi). We also 
assume that 

^r.r _ ^r.r 

where T is some strict open subset of dQ. Write 

Ti = on \ r 

for the inaccessible part of the boundary. The results in this section 
show that in cases where Tj satisfies certain geometric restrictions, it 
is possible to conclude that 

qi = q 2 in tt H (R x O) 

where the sets Oct 2 will be described below. 

Remark. We also obtain results for the conductivity equation by mak- 
ing a standard reduction to the Schrodinger equation. More precisely, 
if lii 72 £ C 2 (Q) are positive functions such that CE' r = C^ 2 ' r , then the 
corresponding DN maps satisfy 

A 7l /|r = A 72 /|r for / G # 1/2 (<9ft) with supp(/) C T. 
Boundary determination |K V84] . |SU88] implies that 
7i|r = 72|r, ^71 | r = d^ 2 \r- 
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Writing qj = Ajj /jj , the relation 

-1/2. / -l/2,s . l-i. 



<->n 



and the above conditions imply that the DN maps A qj for the Schrodinger 
equations satisfy 

A q J\r = A q J\ r for / g tf 1/2 (<9ft) with supp(/) c r. 

Thus C^' 1 " = Cg 2 ' r , and we obtain that 

q x = q 2 in Qn (R x 0). 

Write q = qi = q 2 in (R x O). Then 7^ 2 and 7^ are both solutions 
of (—A + q)u = in fin(lxO) having identical Cauchy data on Y . It 
follows that 71 = 72 in any connected component of Q fl (R x O) whose 
intersection with V contains a nonempty open subset of dVt. 

In the following we will use some general facts on limiting Carleman 
weights from |DKSaU09j . where it was proved that any limiting Car- 
leman weight in R 3 has, up to translation, rotation and scaling, one of 
the following six forms: 

1 II / \ ^ 1 1 \x ~\~ 6]| / if) / • \9\ 

xi, log \x \, a,Tg(x 1 + zx 2 ), 1 — f? , log arg(e (x + iex) ). 

\x\ z \x — ei\ z 

Here 9 G [0, 27r), and the argument function is defined by 

argO) = 2arctan— ^j^— , z G C \ {t G R ; t < 0}. 
\z\ + Re(z) 

It was also proved in [D KSaU09l Section 2] that if if is a limiting 
Carleman weight near (f2,e), then V^y? is a unit parallel vector field 
near (Q, g) where 

-1 IV7 1-2 

g = c e, c = |V e y?| e . 

Furthermore, by the proof of [DKSaU09| Lemma A. 5], if {yi,y') are 
coordinates so that V§y? = d yi and if the coordinates y' parametrize 
a 2-dimensional manifold S such that Vgtp is orthogonal to S with 
respect to the g metric, then the metric has the form 

where g is the metric on S induced by g. 
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3.1. Cylindrical sets. This case corresponds to the limiting Carle- 
man weight ip(x) = X\. Suppose that Q C R x Q , where Qq is a 
bounded domain with smooth boundary in R 2 . Let E be an open 
subset of dfl , and assume that 

r.cMx (an \E). 

If Q has strictly convex boundary, Theorem [2^2] and the result of |Kr09] 
imply that 

q 1 = q 2 in H H (R x O) 

where O is the intersection of Qq with the union of all lines in R 2 that 
have 8Qq \ E on one side. 

The above conclusion holds true also when Qq does not have strictly 
convex boundary. To see this, let Qq CC B CC B where B and B are 
balls. The extensions of the line segments in O to B form a class J 7 
such that any curve in J 7 can be deformed to a point through curves in 
T . It is then enough to extend q% — q 2 by zero to R x B, and to use the 
proof of Theorem 12.21 with M replaced by B together with the result 
[Kr09] . 

3.2. Conical sets. Consider the limiting Carleman weight tp(x) = 
log | x | . Suppose that Q C {£3 > 0}, let (S 2 ,go) be the sphere with 
its standard metric, let S 12 = {u E S 2 ; u 3 > 0}, and let (M ,g ) be a 
compact submanifold of (S+,go) with smooth boundary. Let E be an 
open subset of dM , and assume that 

r,c{ra;r>0,wG dM \ E}. 

We have c = |Vy?| -2 = \x\ 2 and g = \x\~ 2 e, Vgip = x. Choose coordi- 
nates so that 

yi = log|x|, y' = x/\x\. 

The coordinates y' parametrize the manifold S 2 and the metric ^0 on 
S 2 induced by g is just the standard metric g . The discussion in the 
beginning of this section shows that 

Now (M ,g ) is contained in some simple submanifold (M ,g ) of the 
hemisphere (S+,go) (just remove a neighborhood of the equator). Since 
geodesies in S+ are restrictions of great circles, Theorem 12.21 and the 
local injectivity result |Kr09] imply as in Section [37T1 that 



qi = q 2 in Q fl {ruj ; r > 0,u G 0} 
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where O is the union of all great circle segments in such that dM Q \E 
is on one side of the hyperplane containing the great circle segment. 

3.3. Surfaces of revolution. Let OcM 3 \{i;xi < 0}, and consider 
the limiting Carleman weight 

(p(x) = arg(xi + ix 2 ). 

Then V ^=(^,^,0) and 

c = xl + x 2 2 , g = 1 e, Vg(p=(-x 2 ,X!,0). 
x\ + x 2 

We make the change of coordinates valid near Q, 

yi = arg(xi + ix 2 ), y 2 = sjx\ + x 2 2l y 3 = x 3 . 

The coordinates y' parametrize the manifold S = {(xi, 0,x 3 ) ; x\ > 0} 
and Vgtp is orthogonal to S. Furthermore, we may also think of S as 
the set {(0,1/2,1/3) ; V2 > 0}, and the metric on S induced by g is the 
hyperbolic metric go = \e. The discussion in the beginning of this 
section shows that 

«"■•»') = (i w) )■ 

Let (M Q ,g ) be a compact submanifold of S with smooth bound- 
ary, let E be an open subset of <9M . We think of M as lying in 
{(xi, 0,x 3 ) ; x\ > 0}. Now, assume that 

T t c{R e (dM \E);6e(-7r,n)} 

where Rqx = (Rg(xi,x 2 ) t ,x 3 ) t and Rg rotates vectors in R 2 by angle 9 
counterclockwise. That is, we assume that the inaccessible part Tj is 
contained in a surface of revolution obtained by rotating the boundary 
curve <9M \ E. 

Now, the geodesies in S (and, after restriction, also in M ) have 
either the form 

iy2(t),y 3 {t)) = (Rsint, Rcost + a) 

where t G (0, 7r), R > 0, and a G R, or the form {y 2 (t) , y 3 (t)) = (t,a) 
where t > and a G R (these are not unit speed parametrizations). 
In the x coordinates, these are either the half circles in the {x 2 = 0} 
plane given by 

(xi(t) : x 2 (t) : x 3 (t)) = (Rsint, 0, Rcost + a), t G (0, ir) 

or the lines 

(x 1 (t),x 2 (t),x 3 (t)) = (t,0,a), t>0. 
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Enclosing Mq in some ball B in S, the manifold (B, go) is simple and 
Theorem 12.21 and [Kr09] imply as in Section 13.11 that 



where O is the union of all geodesies in S that have dMo \ E on one 
side. 

3.4. Other limiting Carleman weights. So far we have considered 
three of the six possible forms of limiting Carleman weights in M 3 . 
The fourth one, tp(x) = ™, is the Kelvin transform of the linear 
weight, and corresponds to inaccessible parts of the boundary that are 
Kelvin transforms of cylindrical domains. In particular, if part of the 
cylindrical domain is on the hyperplane {x% = 1}, its Kelvin transform 
lies on the sphere centered at (0, 0, 1/2) with radius 1/2, and we recover 
the result of Isakov [Is07j for domains where the inaccessible part is part 
of a sphere. The corresponding results for the remaining two limiting 
Carleman weights do not seem so easy to state and we omit them. 

3.5. Extension of [KSU07] . Let now O C I 3 be a bounded domain 
with smooth boundary, assume that is not in the convex hull of Q, 
and let <p(x) = log \x\. Define 



It was proved in [KSU07j that whenever is a neighborhood of <9Q_U 
dQtan and IV is a neighborhood of dQ + U dQ t an, then 



In particular, and always need to overlap. This result is a 
consequence of the reduction given above for the logarithmic weight, 
Theorem 12.11 (the special case where E = dQ , so that 1^ = 0), and 
injectivity of the ray transform. If dVt tan has zero measure in dQ, then 
Theorem 12.31 allows to improve this result: we have 



In this case, the sets where Dirichlet and Neumann data are measured 
are disjoint, but their union covers all of dQ except for a set of measure 
zero. The result remains true if the roles of dQ + and <9f2_ are changed. 



q x = q 2 in fi n {R e (0) ; 9 G (-vr, vr)} 



dfl± = {xedQ; ±d v <p(x) > 0} 
dVL tan = {i£ dVL ; d v tp(x) = 0}. 
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3.6. Extension of p)7j . The results in [Ts07j stated that C£ r = C T q f 
implies q\ = q 2 in Q if Q C {£3 > 0} and Tj C {£3 = 0}, or if Q C B 
for some ball B and Tj C <9-B. We have already recovered these results 
in Sections 13. II and 1 3.41 since in these cases the local injectivity set O is 
so large that the result q\ = q 2 holds in all of Q. Of course, the results 
above also extend |Is07] since we can conclude at least local uniqueness 
for potentials when the inaccessible part of the boundary satisfies a 
(conformal) flatness condition in only one direction, such as being part 
of a cylindrical set, a conical set, or a surface of revolution. 

We also get global uniqueness if the local injectivity set O is suffi- 
ciently large. For instance, if 

OcRx fi , r.clx (on \ E) 

where f2 is a bounded domain with smooth boundary and E is a 
nonempty open subset of dQ$ , and if the lines in R 2 that have dQo \ E 
on one side cover a dense subset of Qq, then we have q\ = q 2 in Q. One 
example of this situation is if 

fiCtx {(x 2 ,x 3 ) ; x 3 > r)(x 2 )}, Ti C R x {(x 2 ,x 3 ) ! ^3 = ^(^2)} 
where r] : R — > R is a smooth concave function. 

4. Carleman estimate 

Let (M, g) be a CTA manifold, so (M, g) is compact with boundary 
and 

{M,g) CC (R x M ,g), g = c {e®g ). 

Here (Mo, go) is any compact (n — 1) -dimensional manifold with bound- 
ary. We wish to prove a Carleman estimate with boundary terms for 
the conjugated operator e v ^ h (— A g )e~ ip ^ h in M, where ip is the limit- 
ing Carleman weight <p(x) = X\ or (f(x) = —x\, and h > is small. 
Following |KSU07j . it is useful to consider a slightly modified weight 

ip £ = ip + hf £ 

where f £ is a smooth real valued function in M depending on a small 
parameter e, with e independent of h. The convexity of f e will lead to 
improved lower bounds in terms of the L 2 (M) norms of u and hVu. 
On the other hand, the sign of d v (p £ in the boundary term of the Car- 
leman estimate will allow to control functions on different parts of the 
boundary. Of special interest is the set dM tan where d u ip = 0, and in 
this set we have 

d u <Pe\dM tan = hd u f £ . 
We would like to have d v f £ < on dM tan . It is not easy to find a global 
convex function f £ satisfying the last condition for a general set dM tan . 



CALDERON PROBLEM WITH PARTIAL DATA 



17 



However, splitting f e to a convex part whose normal derivative vanishes 
on dM tan and another part which ensures the correct sign on dM tan 
will give the required result. We will use semiclassical conventions in 
the next proof, sec |DKSaU09, Section 4] and [Zwl2j for more details. 
We also write (v,w) = (v,w)l2^, \\v\\ = ||t>|| i 2( M ), and for T C dM 
we write (v,w)r = (v,w) L 2, T y 

Proposition 4.1. Let (M,g) be as above, let <p(x) = ±x\, and let k 
be a smooth real valued function in M so that d v n = —1 on dM. Let 
also q G L°°(M). There are constants e, Cq, h Q > with h < e/2 < 1 
such that for the weight 

ip £ = ip H — f- hn 

e I 

where < h < ho, one has 

h 3 h 2 

— (\d u ip £ \d u u, d v u) dM _ {ips) + — (\\u\\ 2 + || hVu\\ 2 ) 

< \\e^ h (-h 2 A g + h 2 q)(e^ h u)\\ 2 + h 3 (\d^ £ \d„u, d u u) dM+{tpe ) 
for any u G C°°(M) with u\q M = 0. 

Proof. Since (p(x) = ±x± is a limiting Carleman weight in a manifold 
strictly containing M, the computations in the proof of [DKSaU09, 
Theorem 4.1] apply and we can follow that proof. First of all, note 
that 

cV(_A g + q)u = (-A c -i g + g c )(cV u ) 

where q c = cq + c^A^c -1 ^). Thus, by replacing q with another 
potential, we may assume that c = 1 so that g = e © go and <p is a 
distance function in the g metric, i.e. |V ff y| s = 1. 

Let P = -h 2 A g and P , Ve = e^/ h P e~^/ h . Then P 0jlpe = A + iB 
where A and B are the formally self-adjoint operators 

A = -h 2 A g - \Vtp £ \ 2 , B = -2i(V<p £ , hV ■ ) - ihA g ip £ . 

Assume u G C°°(M) and u\qm = 0. We have 

||P ,^|| 2 = {{A + %B)u, {A + iB)u) 

= \\Au\\ 2 + \\Bu\\ 2 + i(Bu, Au) - i{Au, Bu) 

= || Au || 2 + || Bu || 2 + (i[A, B}u, u) - ih 2 (Bu, d u u) dM 

= \\Au\\ 2 + \\Bu\\ 2 + (i[A,B]u,u) - 2h 3 {{d v ^ £ )d v u,d v u) dM . 



Define 



~ / \ hip 2 
<p e (x) = ip + - — . 
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Thus ip £ = (p £ + hn. Let 

A = -h 2 A - |V<^ £ | 2 , B = -2i(V(p e , hV ■ ) - ihAfi £ . 
Since A(p £ = A(p e + hAn and Vy2 £ = V<^ £ + hS/K, we have 

A = A + A e , A e = -h 2 \Wn\ 2 - 2h{V<p £ , V/c), 
B = B + B e , B e = -2ih(Vn, hV ■) -ih 2 An. 
Consequently 

B] = i[A, B] + i[A, B e ] + i[A e , B] + i[A e , B e \. 

Recall from [DKSaU09l p. 143] that 

„ _ Ah 2 ( h \ 2 ~ ~ 
i[A,B] = — \l + -ip\ +hB(3B + h 2 R 

where (3 = (h/e)(l + (h/e)ip)~ 2 and R is a first order semiclassical differ- 
ential operator whose coefficients are uniformly bounded with respect 
to h and e if we assume that h/e < 1/2. Consider now 

i[A, B e ] = i[-h 2 A - |V<^ e | 2 , -2ih{Vn, hV ■) - Hi 2 Ak}. 

It is clear that this equals h 2 Q, where Q is a second order semiclassical 
differential operator whose coefficients are uniformly bounded in h and 
e. The terms z[A e ,-B] and i[A e ,B e ] are better. We thus have 

Ah 2 ( h \ 2 ~ ~ 
i[A, B] = —— ( 1 + —ipj +hB(3B + h 2 Q 

for some Q as described above. It follows that 

Ah 2 ~ ~ 

(i[A, B]u, u) = — 1|(1 + lvp/e)u\\ 2 + h(Bf3Bu, u) + h 2 {Qu, u). 

We will choose ho so small that \h(p/e\ < 1/2 in M for h < h$. Since 
u\bm = 0, integration by parts gives 

~ ~ h 2 ~ 

\h(BfiBu,u)\ < Ci — 1| Bu\\ 2 

Similarly 

\h 2 {Qu,u)\ < C 2 h 2 {\\u\\ 2 + ||Wu|| 2 ). 
Putting this information together, we get 

(i[A,B)u,u) > ^\\u\\ 2 - C^\\Bu\\ 2 - C 2 h 2 {\\u\\ 2 + \\hVuf). 
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Next we revisit the term || Au|| 2 . Let K be a positive constant whose 
value will be specified later. Since u\qm = 0, integration by parts and 
Young's inequality give that 

h 2 \\hVu\\ 2 = h 2 (-h 2 Au,u) = h 2 (Au,u) + h 2 (\V<p £ \ 2 u,u) 
<^\\Au\\ 2 + ^f\\u\\ 2 + C 3 h 2 \\u\\ 2 , 

or 

\\Au\\ 2 > 2Kh 2 \\hVu\\ 2 - K 2 h 4 \\u\\ 2 - 2KC 3 h 2 \\u\\ 2 . 
Also recall that B - B = B e = -2i/i(V/c, hW ■) - ih 2 An. Thus, 
||(5 - B)u\\ 2 < C 4 h 2 (\\u\\ 2 + ||W«|| 2 ). 

Hence 

\\Bu\\ 2 < 2\\Bu\\ 2 + 2C 4 h 2 {\\u\\ 2 + || hVuf) 

and 

||^|| 2 >i||^|| 2 -c 4 /i 2 (|H| 2 + ||w^|| 2 ). 

Putting our estimates together, we obtain 

||^o,^|| 2 > 2Kh 2 \\ hVu\\ 2 - K 2 h A \\u\\ 2 - 2KC 3 h 2 \\u\\ 2 

+ \\\Bu\\ 2 - C 4 h 2 (\\u\\ 2 + \\hVu\\ 2 ) + -IMI 2 - Ci-H^H 2 

2 e e 

- C 2 h 2 {\\u\\ 2 + || hVu\\ 2 ) - 2h 3 ((d u <p £ )d u u, d v u) dM . 
At this point, we choose h so small that 

h 2 1 

e 4 

We also make the choice 



K 



ae 



where a is to be determined. Then for h < h 

\\Po, v M\ 2 > -(M 2 + -\\hVu\\ 2 ) - (C 2 + C 4 )h 2 (\\u\\ 2 + \\hVu\\ 2 ) 
e a 

- -4-H 2 - -— H 2 + hw^f - 2fc s ((^.)^«, d v u) dM . 

e a z e e a 4 



Choose first a = 4C 3 . It follows that 



\\P^M\ 2 > y £ 



1 - 2e{C 2 + C 4 ) 



2^ 

a 2 e 



it 



+ 



a 



- e{C 2 + C 4 ) 



||Wu|| 2 - 2h\{d v <p e )d v u,d v u) dM . 
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Choose next e so that 



mm 



4(C 2 + C 4 )' a{C 2 + C 4 ) 



Finally, choose h so that it satisfies all the restrictions made earlier 
(ho < e/2, h max xeM \(p\ < e/2, and < e/(4Ci)) and additionally 

2^ < 1 
a 2 e ~ 4' 

With these choices, we have 

h 2 h 2 _ 
||-Po, Ve m I| 2 > ^-|M| 2 H II^Vm|| 2 - 2h 3 ((d u (p £ )d u u, d v u) 9M - 

o£ OLE 

Adding a potential gives 

H^V^II 2 < 2||(P 0lVe + h 2 q)u\\ 2 + 2/ i 4 ||g||i 0O(M) || M || 2 . 

Choosing an even smaller value of ho depending on ||<?||l°°(m) if neces- 
sary, we obtain for < h < h that 

IK^V + h 2 q)u\\ 2 > ^(\\u\\ 2 + || hVuf) - 2h 3 ((d^ £ )d„u, cUW 

Finally, we replace u by e ip2 ^ 2£+K u, where u G C°°(Q) and u\qq = 0, 
and use the fact that 

1/C < e^ ,2e+K < C, |V(e^ 2/2e+K )| < C on M. 

The required estimate follows. □ 

We now pass from ip e to (p in the boundary terms of the previous 
result, making use of the special properties of tp e on DM. Note that 
the factor h A in the boundary term on {x G dM ; —5 < d u ip(x) < h/3} 
below is weaker than the factor h 3 in the other boundary terms. This 
follows from the fact that d u (p £ = hd v n = —h in the set where d u ip 
vanishes, so one only has the weak lower bound. 

Proposition 4.2. Let (M,g) be as above, let q G L°°(M), and let 
(p(x) = ±Xi. There exist constants Co, ho > such that whenever 
< h < ho and 5 > 0, one has 

Sh 3 /i 4 h 2 

7^ll^lli 2 ( { a^<^}) + 7rll^ll! 2( {-^^<V3}) + 7r(ll«ll 2 + II^VM|| 2 ) 

Uo Oo Oo 

< \\e^ h (-h 2 A g + h 2 q)(e-^ h u)\\ 2 + h 3 \\d u u\\h i{ 8^> h/ 3 }) 
for any u G C°°(M) with u\qm = 0. 
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Proof. Note that 

d v ip e = ^1 + ~(pj d v ip + hd v n = ^1 + jv 9 ^ d v ip - h. 

We choose ho so small that whenever h < ho, one has for x E M 

1 , h , , 3 
- < 1 + ~<z>(z) < -. 

On the set where d u (p(x) < —5, we have 

|^ £ | > 5/2. 
If —5 < d v ip < h/3, we use the estimate 

\d v (p e \ > h/2. 

Moreover, \d u (p e \ < C on DM. Since {d v ip < h/3} C {d u <p £ < 0} and 
{d u <p e > 0} C {dyip > h/3}, the result follows from Proposition 14.11 
after replacing Co by some larger constant. □ 

We can now obtain a solvability result from the previous Carleman 
estimate in a standard way by duality (see |BU01[ IKSU07t INS10] ). 
There is a slight technical complication since the solution will be in L 2 
but not in H 1 . To remedy this we will work with the space 

H Ag (M) = {uE L 2 (M) ; A g u E L 2 {M)} 

with norm || m || h a = IM|l 2 + ||Aw||l2. As in jBUOlj . we see that Ha(M) 
is a Hilbert space having C°°(M) as a dense subset, and there is a 
well defined bounded trace operator from H A (M) to H~ l l 2 (dM) and a 
normal derivative operator from H A (M) to H~ 3 / 2 (dM). We also recall 
that 'due H A (M) and u\ dM E H*/ 2 (dM), then u E H 2 (M). 

Proposition 4.3. Let (M,g) be as above, let q E L°°(M), and let 
<f(x) = ±x±. There exist constants Co, To > such that when t > tq 
and 5 > 0, then for any f E L 2 (M) and /_ E L 2 (S- U So) there exists 
u E L 2 (M) satisfying e TLp u E H Ag (M) and e TLp u\ dM e L 2 (dM) such 
that 

e-^(-A g + q){e T *u) = f in M, e^u\ s _ uSo = e^/_, 

and 

\\u\\ lH m) < C (r- l \\f\\ LHM) + (Sry^Wf-lsJ^ + ||/_| So |U 2 ( So) ). 
Here S± and So are the following subsets of dM: 
S- = {d uV < -5}, So = {-5 < dy V < l/(3r)}, S + = {d u ip > l/(3r)}. 
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Proof. Write Lv = e Tip (-A g + q)( e - Tlp v) and r = 1/h, r = l/h . We 
rewrite the Carleman estimate of Proposition 14.21 as 

(5r) 1/2 \\d u v\\ L 2 {s _) + \\d v v ||l2(5 ) + t\\v || + ||Vv|| 

<C7 ||Li;||+C7or 1/2 ||^|| L2(5+) . 

This is valid for any 5 > 0, provided that t > tq and t> G C°°(M) with 
f |aM = 0. 

Consider the following subspace of L 2 (M) x L 2 (S + ): 

X = {(Lv, d u v\ s+ ) ; v G L7°°(M), = 0}. 

Any element of X is uniquely represented as (Lv, d v v\g + ) where v\qm = 
by the Carleman estimate. Define a linear functional I : X — > C by 

l(Lv,d v v\ s+ ) = (vJ) L 2 {M) - (d u vJ„) L 2 {s _ uSo) . 

By the Carleman estimate, we have 

\l(Lv,d v v\ S+ )\ < \\v\\\\f\\ + ||9 v u||La(s_)||/_||La(s_) 

+ \\dvv\\ L 2 {So) \\f_\\ L 2 iSo) 

< Co^WfW + (ST)-^\\f-\\^8-) + ||/-||L»(flb)) 

x (ii^n + r 1 / 2 !!^!!^^)). 

The Hahn-Banach theorem implies that I extends to a continuous linear 
functional J: L 2 (M) x T~ 1 f 2 L 2 (S + ) ->■ C such that 

IKII < C (T- l \\f\\ + (^r)- 1 / 2 !!/^!^^) + ||/_|U 2( 5 0) ). 

By the Riesz representation theorem, there exist functions u G L 2 (M) 
and u + G L 2 (S+) satisfying l(w,w + ) = (w,u) L 2 {M) + (w + ,u + ) L 2 {s+) . 
Moreover, 

IMU 2 (Af) + t ~ 1/2 \\u+\\l 2 {S+) 

< C (r- l \\f\\ + (Srr^Wf^s^ + Wf-hHSo))- 
live C°°(M) and v\ aM = 0, we have 
(4.1) (Lv,u) L 2 (M) + (<9^,m+) L 2 (5+) = (v,f) L 2 {M) 

- (d v v,f-) L 2 (s _ uSo) . 

Choosing v compactly supported in M mt , it follows that L*u = f, or 

e~ TLp (-A g + q)(e Tlp u) = f in M. 

Furthermore, e Tip u G H&(M). 

If w,v G C°°(M) with v|aM = 0, an integration by parts gives 

(Lv,w) = -(e~ Ttp d v v,e Ttp w) L 2^ M) + (v,L*w). 
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Given our solution u, we choose Uj G C°°(M) so that e Ttp Uj — > e Tip u 
in H&(M). Applying the above formula with w = Uj and taking the 
limit, we see that 

(Lv,u) = -{e-^duV, e T ^u) L 2 (gM) + (v,L*u) 

for v G C°°(M) with v\qm = 0. Combining this with (14. ip . using that 
L*u = f, gives 

(d v v , /_)£3(5_us ) + (^,«+)l 2 (5+) = (e" Tip <9^,e Tl/, u) L 2 (aM) . 
Since d u v can be chosen arbitrarily, it follows that 

e Tlp u\s„us = e Tip f-, e Tip u\ s+ = e TLp u+. 
We also see that e T *u\ dM G L 2 (dM). □ 

5. Reflection approach 

In the previous section, we employed Carleman estimates and dual- 
ity to obtain a solvability result (Proposition 14. 3[) that will be used to 
produce correction terms in complex geometrical optics solutions with 
prescribed behavior on parts of the boundary. In this section we give 
an alternative approach to the construction of correction terms van- 
ishing on parts of the boundary. The method is based on a reflection 
argument. We extend the method of [Is07j, which dealt with inacces- 
sible parts that are part of a hyperplane, to the case of inaccessible 
parts that are part of the graph of a function independent of one of 
the variables. The results are less general than the ones in Section HI 
and for simplicity will only be stated for domains in IR 3 with Euclidean 
metric, but on the other hand the method is constructive and is based 
on direct Fourier arguments in the spirit of [KSUlla] , |KSUllb"] . 

Let fl C 8 3 be a bounded domain with smooth boundary, and as- 
sume that 

n C R x {(x 2 ,x 3 ) ; x 3 > 77(x 2 )} 

where r\ : K — > M is a smooth function. Also assume that Tq is a closed 
subset of dfl such that 

T C R x {(x 2 ,x 3 ) ; x 3 = i](x 2 )}. 

We will show that if one has access to suitable amplitudes of complex 
geometrical optics solutions that vanish on Tq, then it is possible to 
produce correction terms that also vanish on r . 

Proposition 5.1. Let Q and Tq be as above, and let q G L°°(Q). 
There are C , r > such that for any r with |r| > r and for any 
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m E H 2 (fl) with m|r = 0, the equation (—A + q)u = in Q has a 
solution u G H 2 {Q) of the form 

u = e~ TXl (m + r) 

such that r|r = and 

\\r\\i?<ss> < %e^(-A + q)(e-^m)\\ L 2 (u) . 

The proof involves a reflection argument that reduces the construc- 
tion of the correction term to the problem of solving a conjugated 
equation with anisotropic metric, 

e TXl (-Ag + q) {e~ TXl r) = f in R x tt 

where Cl C R 2 is a bounded open set and g is a metric of the form 

and where go is smooth for 1/3 ^ but only Lipschitz continuous across 
{y 3 = 0}. In three and higher dimensions, it is not known how to handle 
equations of this type with general Lipschitz coefficients in the second 
order part (the case of C 1 coefficients, and also Lipschitz coefficients 
with a smallness condition, is considered in [HTllj ). However, in our 
case the singularity of g only appears in the lower right block g , and 
this turns out not to be a problem. 

The following is an analogue of [KSUlla] Proposition 4.1], the main 
difference being that the transversal metric is only Lipschitz. (With 
correct definitions, one could easily deal with L°° transversal metrics 
as well, but then the solution would only be in H^_ S (T).) Here we write 
(xi,x') for coordinates in T = R x M , and for (Jetwe consider the 
spaces 

H/|Ul(T) = \\(xi) S f\\L*(T), \\f\\HUT) = ll/IUi(T) + IM/IU i( T) 

with (t) = (l+t 2 ) 1/2 , and similarly for Hj(T). We also write Spec(-A 3o ) 
for the set of Dirichlet eigenvalues of the Laplace-Beltrami operator 
-A go in (M ,g ). 

Proposition 5.2. Let T = R x M with metric g = e © go, where 
(Mo, go) is a compact oriented manifold with smooth boundary and g 
is a Lipschitz continuous Riemannian metric on M Q . Given any q e 
L^ mp (T) and any 5 > 1/2, there are constants C ,r > such that 
whenever 

|t| > t and t 2 £ Spec(— A 50 ), 
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the equation 

e TXl (-A g + q)(e- TXl r) = f inT 

has a unique solution r G H\ & (T) with t\qt = for any f G L 2 (T). 
Moreover, r G H^_ S (T), and one has the bounds 

Co 

\\r\\i? t p>) < rr II/IUk t )' H r lltfI 4 Cr) < Coll/llz-JCT)- 

Proof. The proof is almost exactly the same as the proof of |KSUlla"t 
Proposition 4.1], and we only give the main idea. Since A 9 = d 2 + A 9o , 
the equation that we need to solve is 

(-d 2 xi + 2rd Xl - A go -r 2 + q)r = f in T. 

It is enough to consider q = 0. The standard argument based on weak 
solutions shows that even when go has very little regularity, there is an 
orthonormal basis of L 2 (M ) consisting of Dirichlet eigenfunctions of 
-A ff0 , 

-A go (pi = Xicpi in M , (pi G #o( M o), 

where < Ai < A2 < A3 <...—»■ 00 are the Dirichlet eigenvalues of 
-A ff0 in M . 

Considering the partial Fourier expansions 

00 00 

r(x u x') = ^ O^iOA f(xi,x') = ^2 K x u l)<Pift), 
1=1 1=1 

it is enough to solve 

{-d 2 Xi + 2rd Xl +\- r 2 )r{ -,/) = /(.,/) in R for all I. 

The condition r 2 ^ Spec(— A go ) allows to solve these ordinary differ- 
ential equations by the Fourier transform as in |KSUlla"t Section 4], 
and the estimates given there imply that one obtains a unique solution 
r G H^_ S (T) with r|^T = satisfying the required bounds. Elliptic H 2 
regularity also works with Lipschitz go, and the argument in [ KSUlla 



Section 4] gives that r G H 2 S (T). □ 

Proof of Proposition \5.1\ We begin by flattening T via the map 

$ : R 3 -»■ M 3 , (x 1 ,x 2 ,x 3 ) H> (a*, 0:2,2:3 -17(0:2)). 

Let f2 = $(f2), write y for coordinates in Q, and let i? be the reflection 

R(yi,y2,ys) = (2/1,2/2, —2/3)- 

Note that Q C {j/3 > 0}. Consider the reflected domain f2* = R(Q), so 
ft* C {2/3 < 0}, and let U the double domain Cl U $(r ) int U Q*. 



26 



CARLOS KENIG AND MIKKO SALO 



Let \1/ = let g = \P*e be the metric in Q that is the pullback 
of the Euclidean metric in f2, let q = ^*q, and let rh = ty*m. In the 
double domain U, we use even reflection to define the quantities 

9, V3 > 0, A \ q, y 3 > 0, 



* R*~g, y 3 <0, q { R*q. !,,<(). 
and odd reflection to define the amplitude 



m 



1 m, y 3 > 0, 



2 



-\R*m, y 3 < 0. 
Since the flattening map <3> leaves X\ intact, we have 



g(yi,y') 



l o 

o g (y') 



where go is a Lipschitz continuous metric only depending on y 2 and y 3 . 
(In fact, g and go are well defined in {y 3 > 0} by the flattening map 
$ and the Euclidean metric in {x 3 > 77(^2)}. ) Also, q G L°°(U), and 
rh G H 2 (U) by the boundary condition m|r = and by the properties 
of odd reflection. 

We wish to find f G H 1 ^) satisfying 

e™ 1 (-A« + ?)(e- TXl f) = / 
where / = — e TXl {— A§ + q){e~ TXl m). Now 

ll/IU^) = \\f\y { a) + WfWmn*) 

= ||^(e"' 1 (-A + ?)(e-™ 1 m)|| £a(fl) 

+ || J R*^(e^(-A + g)(e-^m)|| L2( ^ ) 
<C7||e™ 1 (-A + 9)(e-™ 1 m)|| L2(n) . 

Choose a bounded open set &o C R 2 such that 

U CC K. x Q , 

and let go be the metric in Qo that is the even extension of go from 
{y 3 > 0} to Cl . Then g is smooth for y 3 7^ and Lipschitz continuous 
across {y 3 = 0}. Extending g to 1 x fi using the block structure and 
extending q and / by zero to R x Q , it is enough to find a solution 
r G Hf oc (R x Ct ) of the equation 

(5.1) e rai (-A§ + q)(e~ TXl r) = f inRx Q . 
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Such a solution may be found by Proposition 15.2} and denoting by f 
its restriction to U we have 

C 

\\r\\i*(u) < t-tII/IU 2 ^)- 

Define now 

u = e~ TXl (rh + f) in U 

and 

u = u — R*u in Q. 

Then (— A^ + q)u = in U, and (— A^ + q)u = in Q by the definition 
of g and q and using that u G H 2 {U). We also have 

u = e"™ 1 (m - R*m + r - R*r) in fi. 

But here rh — R*rh\^ = rh by the definition of rh. Consequently, if we 
define u = then (—A + q)u = in Q and 

u = e~ TXl (m + r) in Q, 

where r = $*(f — i?*f) satisfies 

C 

lkl|jy»(n) < C||^IIl2(c/) < o Wfh^u) 

<^\\e^(-A + q)(e-^m)\\ L2m . 

This proves the result. □ 

Note how the odd reflection of the amplitude m in the proof ensured 
that the solution obtained by reflection is not the zero solution. We 
also remark that under certain conditions, the arguments in Sections [6] 
and [7] allow to construct amplitudes m vanishing on a part r as above. 

6. Local uniqueness on simple manifolds 

In this section we prove Theorems I2.1H2.31 In these results the 
transversal manifold is assumed to be simple and we only use non- 
reflected geodesies. This case already illustrates the main features of 
the approach, and we can use a quasimode construction that is much 
easier than the Gaussian beam one used for non-simple transversal 
manifolds and reflected geodesies. 

The first observation is the usual integral identity. 

Proposition 6.1. IfV D ,V N C dM are open and if ' C v g °f N = C v g °f N , 
then 




U1U2 dV g = 
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for any Uj G H& g (M) satisfying (— A g + qj)uj = in M and 

supp0i|,9 A /) C T D , supp(u 2 \ dM ) C IV 

Proof. Let Uj be as stated. Since C^®f N = C^°f N , there is a function 
u 2 G Ha(M) with (—A + g 2 )w 2 = in M, supp({£ 2 |sAf) C T D) and 

{ui\v D ,d v ui\T N ) = {u 2 \r D ,d u u 2 \r N ). 
Using that u\, u 2 and u 2 are solutions, we have 

/ (qi - Q2)u\u 2 dV = / [(Aux)u 2 - u 1 (Au 2 )) dV 



[(A(mi - u 2 ))u 2 - (ui - u 2 ){Au 2 )\ dV. 

' M 

Now Ui — u 2 \qm = 0, so in fact u\—u 2 G H 2 (M) by the properties of the 
space Ha(M). Recall also that C°°(M) is dense in H&(M) and that 
u 2 \ 9 m e H- 1/2 (dM) and d u u 2 \ dM G H- 3/2 (dM). These facts make it 
possible to integrate by parts, and we obtain that 

/ (?i - <l2)uiu 2 dV = [(d v (ui - u 2 ))u 2 - Oi - u 2 )(d u u 2 )} dS 
Jm JdM 

in the weak sense. The last expression vanishes since d v (u\— u 2 )\r N = 

and supp(M 2 |9A/) C T n . □ 

The next result will be used to pass from the metric g = c(eQ)go) to 
the slightly simpler metric g = e © go. 

Lemma 6.2. Let c be a smooth positive function in M. Then u G 
Ha 9 (M) satisfies (— A g + q)u = in M if and only if u G H&.(M) 
satisfies (—A g + q)u = in M, where 



n — 2 . n — 2 , n — 2 , 



g = c g, u = c 4 u, q = c(q — c 4 A g (c 4 )). 
Proof. This follows from the identity for v G C°°(M), 

n+2 , n—2 n-2 n-2 

c 4 (-A g + q)(c 4 v) = (-A c -i g + c(q-c 4 A g (c 4 )))v, 
upon approximating u or u by smooth functions. □ 

Proof of Theorem ] 2. 1[ Let g = e © g and q~j = c(qj — A g (c~ n ~^)). 
Let A be a fixed real number, and consider the complex frequency 

s = r + i\ 

where r > will be large. We look for solutions 

ui = e- SXl (v s (x')+ ri ) } 
u 2 = e SXl {v s {x')+r 2 ), 
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of the equations (— A^ + qx)U\ = 0, (— A g + ^2)^2 = in M. Here 
v s will be a quasimode for the Laplacian in (Mo, go) that concentrates 
near the given geodesic 7. Next we will construct a suitable solution 
Ui, and the case of u 2 will be analogous. 

Since A g = df + A go , the function ui is a solution if and only if 

(6.1) e s ^(-A g + giXe-^n) = -(-A 9o + & - s 2 )v s (x') in M. 
We want to choose v s G C°°(M ) to satisfy 

(6.2) \\v s \\ LHMo) = 0(1), ||(-A go - s 2 )v s \\ L2[Mo) = 0(1) 
as r — > 00. Looking for f s in the form 

w s = e^a 

where ^,o6 (^(Mo), a direct computation shows that 

(~ A 90 - 

= e*** (s 2 [|#|p - 1] a - is [2<#, d ■ ) So + A, ^] a - A go a) . 

Since (Mo, go) is simple, it is easy to find tp and a so that the expres- 
sions in brackets will vanish and that the resulting quasimode v s will 
concentrate near the geodesic 7. To do this, let (Mo, go) be a simple 
manifold that is slightly larger than (M , go), extend 7 as a geodesic in 
M , and choose e > such that 7|(-2e,o)u(L,L+2 £ ) stays in M \M (this 
is possible since 7 is nontangential). Let u = j(—e) G M \ M , and 
let (r,6) be polar normal coordinates in (Mo, go) with center u. Then 
7 corresponds to the curve r 1— > (r,9o) for some fixed 8q G S^ -2 . We 
will choose 

i/)(r, 9) = r, 
a(r,6) = \go(r,9)\- l '%(6) 

where \go\ is the determinant of go, and b is a fixed function in C°°(S n ~ 2 ) 
that is supported so close to #0 such that v s \qm \e = 0. With these 
choices, we have as in |DKSaU09] 

(_ Ago - s 2 )v s = -e^A go a. 

Thus v s satisfies the estimates (16.21) . and also the estimate 

||fs|U°°(Mo) = 0(1). 

We now go back to (16. ip . and look for a solution in the form 77 = 
e i\x lr ^ w here satisfies 

(6.3) e TX1 (-A g + qi)(e- TXl r[) = f in M 
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with 

f = -e- i ^(-A go + q 1 -s 2 )v s (x'). 

We also want to arrange that supp^i]^) C Tc where Yd 3 <9M_ur a . 
For this purpose, let S > be a small number to be fixed later, let S± 
and So be the sets in Proposition 14.31 with Carleman weight (p(x) = 
—Xi, define 

V s = {x e S_ U S ; distdA/O, r<) < 5 or x e dM + }, 
T 5 a = (S_ U So) \ 
and impose the boundary condition 

(6.4) e r *r[\s_uSo=e T ' p f- 
where 

f -e- iXxi v s (x'), on 
" \ 0, on rt 

Note that dM + U dM tan (these sets refer to the weight X\) is in the 
interior of S_ U So in DM. 

We have seen that ||/||l 2 (A:0 = 0{1) as r — > oo. We also have 

f-\ dMtan — 0) 

since /_| r * naMtan = by definition and f-\dM tan nv s = for sufficiently 
small 5 > by the construction of v s and using that ^ C Rx (dM \E). 
Since ||/_||x-(5_u5 ) < 1, we have 

Wf-\\V(S-) < '{{fa! > 5}) 

and 

ll/-ll^(So) £ ^({-V(3r) < d v x x < 0} U {0 < d v x x < 5}), 

where a is the surface measure on DM. It follows from Proposition 
14.31 that the equation (16. 3p has a solution r[ satisfying the boundary 
condition (16.41) . and having the estimate 

\K\\l H m) < r- 1 + (Srr^add^ > 5}) 

+ <x({-l/(3r) < d v x x < 0}) + a({0 < d v x x < 5}). 

The implied constants in the previous inequality are independent of r 
and S. By the basic properties of measures, for some constant C > 
we have 

\K\\l H m) < C [t- 1 + (St)- 1 ' 2 + 0^(1) + o^o(l)] • 
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Given e > 0, we first choose S so that CqOs->o(1) < e/2. After this, we 
choose r > so large that C (r- 1 + (c^r)" 1 / 2 + o^^l)) < This 
shows that 

lim |K(-;r)|| L 2 (M) = 0. 

Choosing r' x as described above and choosing rj = e^ 1 ^, we have 
produced a solution e H&~(M) of the equation (— + q~i)ui = in 
M, having the form 

mi = e- SXl (v s (x') +n) 

and satisfying 

supp(£ti| aA/ ) C T D 

and || ri ||i,2(m) — o(l) as r — > oo. Repeating this construction for the 
Carleman weight <f(x) — x%, we obtain a solution u 2 € H^ S (M) of the 
equation (— + ^2)^2 — in M, having the form 

u 2 = e SXl {v s {x')+r 2 ) 

and satisfying 

supp(m 2 |om) C T n 
and ||?~2||l 2 (a/) — o(l) as r — >• 00. 

n _2 

Writing = c Lemma 16.21 shows that Uj G H^ g (M) are 

solutions of (— A g + (71)1*1 = and (— A g + q2)u 2 = in M. Then 
Proposition 16.11 implies that 

(qt - q 2 )u 1 v^dV g = 0. 

M 

We extend q\ — q 2 by zero to R x Mq. Inserting the expressions for Uj, 
and using that dV g = c n l 2 dx\ dV go (x'), we obtain 

/ / (qi-q2)ce- 2iXxi (\v s (x')\ 2 + v s r 2 - + - s r 1 +r 1 rl)dx 1 dV go (x') = 0. 

./Mb ./-oo 

Since ||rj||i2( M ) = o(l) as r — > 00 and since <iVg = |^o| 1 ^ 2 ^'"^ m the 
(r, 0) coordinates, it follows that 

e~ 2Xr {c{qi - <? 2 )r(2A, r, 0) |6(6>) | 2 dr ^ = 0. 

s™- 2 ./o 

Varying 6 in C°°(S n ~ 2 ) so that the support of b is very close to 9 , this 
implies that 

00 

e- 2Xr (c( qi -q 2 )n2\,r,9 )dr = 0. 
Since 7 was the curve r h-> (r, 9), this shows the result. □ 
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Proof of Theorem \2.2i Suppose that the local ray transform is injective 
on O and O D dM C E. By Theorem l2.lt we know that 

(6.5) [ L e- 2Xt (c( gi - g 2 )R2A, j{t)) dt = 

Jo 

for any nontangential geodesic 7 in O. Setting A = and using local 
injectivity of the ray transform, we obtain that 

(c(<Zi - <Z 2 )R0, -) = 0inO. 

Going back to (16. 5p and differentiating this identity with respect to A, 
and then setting A = and using the vanishing of (c(gi — q 2 ))"(0, • ) on 
O, it follows that 

L ^[(c(gi-g 2 ))1 (0, 7 (t))<ft = 0inO 

for any nontangential geodesic in O. Local uniqueness for the ray 
transform again implies that 

|f [(c(fc - (0, -) = 0inO. 

Iterating this argument by taking higher order derivatives of (16.51) 
shows that 

J^) -<&))! (0, -) = 0inO 

for any fc. Since c(qi — q 2 ) is compactly supported in xi, its Fourier 
transform is analytic and we have 

(c(?i - 92)) A (A, • ) = in O for all A G R. 

Inverting the Fourier transform and using that c is positive, we obtain 
that q 1 = q 2 in M n (R x O). □ 

Proof of Theorem \2.3[ Since (M,g) is admissible, we may assume that 

(M,<7) C (lxM 0)? ), & = c(e©^ ) 

where (M ,g ) is simple. The argument is very similar to the proof 
of Theorem 12.14 an d we only indicate the required changes. Up to 
the formula (16.31) . the only change is that there is no restriction on 
b G C°°(S n ~ 2 ) (we do not require v s to vanish on any part of the 
boundary). The function r[ is obtained as a solution of (16.31) . but this 
time we want that supp(-ui|aM) C <9M_. Fix 5 > 0. The boundary 
condition for u\ is (16. 4p . where /_ is chosen to be 

/_ = -e~ iXxi v s (x') onS_US . 
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We use Proposition 14.31 to solve for r[. We have ||/||l2( M ) = 0(1), 
and the bound < 1 implies 

\\f-\\»(S-) S > 8}) 

and 

II/-IIl»(s ) £ ^({"V(3r) < d v x x < 0})+a(dM tan )+a({0 < d v x x < 5}). 
Now we use that 

a(dM tan ) = 0. 

This shows that we obtain the same estimate for r[ as before: 
KWlhm) < Co [r- 1 + (St)- 1 / 2 + 0^(1) + o 5 _> (l)] • 
We can now continue as in the proof of Theorem 12.11 to conclude that 

/V 2A *(c( gi -g 2 )r(2A,7(t))^ = 
Jo 

for any A G K and for any nontangential geodesic in (Mo, go). The 
geodesic ray transform (with zero attenuation) is injective in (Mo, go) 
[Sh.94] . Following the proof of Theorem 12.21 but now using all the 
nontangential geodesies in (Mo, go), shows that q\ = q<i in M. □ 

7. QUASIMODES CONCENTRATING NEAR BROKEN RAYS 

In this section, to simplify notation, we write (M, g) instead of 
(Mo, go) and we assume that (M,g) is a compact oriented Riemannian 
manifold having smooth boundary and dim(M) = m > 2. Suppose 
that E is a nonempty open subset of DM, and let R = DM \ E. We 
think of E as the observation set where geodesies can enter and exit, 
and R is the reflecting set. In the Calderon problem with partial data 
we are led to consider attenuated broken ray transforms, where one 
integrates a function on M over broken geodesic rays that enter M at 
some point of E, reflect nontangentially at points of R, and then exit 
M at some point of E. The reflections will obey the law of geometric 
optics, so that a geodesic hitting the boundary in direction v will be 
continued by the geodesic in the reflected direction v — v — 2(v, v)v. 

Given a slightly complex frequency s = r + iX, we will construct 
corresponding quasimodes, or approximate eigenf unctions, that con- 
centrate near a fixed nontangential broken ray. 

Proposition 7.1. Let 7 : [0, L\ — > M be a nontangential broken ray 
with endpoints on E, and let A be a fixed real number. For any K > 
there is a family {v s ; s = r+i\, r > 1} in C°°(M) such that as t — >■ 00 

\\(-A g - s 2 )v s \\ L 2 (M) = 0(t~ k ), \\v s \\ L 2 {m) = 0(1), 
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the boundary values of v s satisfy 

IHU 2 (i?) = o(t~ k ), \\v s \\ L 2 idM) = 0(1), 

and for any ip G C(M) 

/ \v T +i\\ 2 ip dV g ->■ / e~ 2Xt tp(j(t)) dt as r ->■ oo. 
Jm Jo 

Let us begin by proving this result in the special case -E = (9M, so 
that R = and one does not need to worry about reflected rays. The 
next three preparatory lemmas describe a modified Fermi coordinate 
system that is very useful in this construction. 

Lemma 7.2. Let (M,g) be a compact manifold without boundary, and 
let 7 : (a,b) — > M be a unit speed geodesic segment that has no loops. 
There are only finitely many times t G (a, b) such that 7 intersects itself 
at 7(f). 

Proof. Since 7 has no loops, (7(£),7(£)) = (7(f), 7(f)) implies t = t! . 
The first observation is that 7 can only self-intersect transversally, since 
also (7(t),7(t)) = (j(t'), —j(t')) implies t — t' (if this would happen 
for t < t', then by uniqueness of geodesies 7(^) = — t(^-) which is 
impossible). Next note that if r is smaller than the injectivity radius 
of (M,g), then any two geodesic segments of length < r can intersect 
transversally in at most one point (locally geodesies are close to straight 
lines). Partitioning (a, b) in disjoint intervals { Jk}k=i °f length < r, we 
have an inject ive map 

{(t, t') G (a, bf ; t < t' and 7(f) = 7^')} 

l— ^ {(^,1) G {1, . . . , K} 2 ; t G Jk,t' G J{\. 

Consequently, 7 can only self-intersect finitely many times. □ 

Lemma 7.3. Let F be a C 1 map from a neighborhood of (a, b) x {0} in 
]R n into a smooth manifold such that -F|( a ,6)x{o} is injective and DF(t, 0) 
is invertible for t G (a, b). If[a , b ] is a closed subinterval of (a, b), then 
F is a C 1 diffeomorphism in some neighborhood of [ao, bo] x {0} in M. n . 

Proof. For any t G [a ,&o]> the inverse function theorem implies that 
there is e t > such that F\(t-3 £tt t+3e t )xB 3s (0) is a C 1 diffeomorphism. 
Since [ao, &o] is covered by the intervals (t — e t , t + e t ), by compactness 
we have [a ,b ] C l>f =1 (tj - e j: tj + Sj) where F|( tj ._ 3ej . )tj . +3ej .) >< B3 ej .(o) is 
bijective. We can further assume (upon throwing away or shrinking 
some intervals if necessary) that the intervals Ij = (tj — £j,tj + Sj) 
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satisfy Ij fl Ik — unless \j — k\ < 1. Since 7(t) = F(t, 0) is injective, 
we also have 'y(Ij) fl ^(h) = unless |j — k\ < 1. 

Fix a Riemannian metric in the target manifold, and define 

5 = inf {dist( 7 (7 j ), 7 (7 fe )) ; \j -k\>2}>0. 

Let Uj = Ij x B £ (0), where e < min{ei, . . . , en} is chosen so small that 
F(Uj) C {q; dist(g,7(7 i )) < 5/2}. Then F(Uj) H F(E/*) = unless 
|j — k\ < 1. Define 

U = U? =1 Ui. 

To show that F|j/ is a C 1 diffeomorphism, it is enough to check in- 
jectivity. If F(t,y) = F(t',y') for (t,y), (t',y') e C/, then necessarily 
(t, y) G C/j, (t',y f ) G C/fe where |j— A;| < 1. We may assume that £j > e^. 
Since F\ {t ._ 3£j:t . +3£ . )xB3e . {0) is bijective, we obtain (t,y) = (t',y'). □ 

Lemma 7.4. Lei (M, (?) be a compact manifold without boundary, and 
assume that 7 : (a, b) — > M is a unit speed geodesic segment with no 
loops. Given a closed subinterval [ao,&o] of(a,b) such that 7|[ ao ,&o] se tf~ 
intersects only at times tj with a < t\ < . . . < t N < b (sett = a and 
tN+i = °o), there is an open cover {(Uj, (fj)}^^ 1 o/ 7 ([a 0) &o]) consisting 
of coordinate neighborhoods having the following properties: 

(1) ipj(Uj) = Ij x B where Ij are open intervals and B = B(0, 5) is 
an open ball in M n_1 where 5 can be taken arbitrarily small, 

(2) <p j ('y(t)) = (t,0)fort el t 

(3) tj only belongs to Ij and Ij fl Ik — unless \j — k\ < 1, 

(4) (pj = if k on (fj^ilj D h) x B). 

Further, if S is a hypersurface through 7(a ) that is transversal to ■j(a ), 
one can arrange that the map y h-> ( / 9Q 1 (a ,y) parametrizes S near 
7(a ). 

Proof. We will use modified Fermi coordinates, constructed as follows. 
Let {i>i, . . . , v n -i} be an orthonormal set of vectors in T 7 ( ao )M such that 
{j(ao),vi, . . . , i> n -i} is a basis. (The case where {j(ao),vi, . . . , f n -i} 
is an orthonormal basis corresponds to the usual Fermi coordinates.) 
Let E a (t) be the parallel transport of v a along the geodesic 7. Since 
i(t) is also parallel along 7, the set {i(t), Ei(t), . . . , _E n _i(t)} is a basis 
of T 7(t) M for t G (a, 6). 
Define the function 



F : (a, 6) x R"- 1 ^ M, F(i,y) = exp 7(t) (y Q F a (t)). 
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Here exp is the exponential map in (M, g) and a, [5 run from 1 to n — 1. 
Then F(t, 0) = j(t) and (with e a the ath coordinate vector) 

—F(t,se a )\ s=0 = E a (t), ^F(t,0) = j(t). 

Thus F is a C°° map near (a, 6) x {0} such that DF(t, 0) is invertible 
for t e (a, 6). 

In the case where 7 does not self- intersect, -F|( afe ) x { } is injective and 
Lemma 17751 implies the existence of a single coordinate neighborhood of 
7([a ,&o]) so that (1) and (2) are satisfied (then (3) and (4) are void). 
In the general case, by Lemma 17721 the geodesic segment 7|[ ao ,M oru y 
self- intersects at finitely many times tj with a < t\ < . . . < < b . 
For some sufficiently small 5, 7 is injective on the intervals (a,ti — 5), 
(ti — 2(5, t 2 — 5), . . . , (t n — 25, b) and each interval intersects at most two 
of the others. Restricting the map F above to suitable neighborhoods 
corresponding to these intervals (or slightly smaller ones) and using 
Lemma 17731 we obtain the required coordinate charts with tpj = F~ l \u r 

Let 5* be a hypersurface transversal to 7(00), and choose some para- 
metrization y 1— > q(y) of S near 7(a ) satisfying -^q(se a ) = v a . We will 
form a new chart (Uq, y5o) by modifying (Uq, y?o) so that y i-)- (p^ 1 (ao, y) 
parametrizes S near 7(00). 

We may assume that a = 0, and write F = ip^ , F = (p^ . It is 
enough to choose F = F Q o $, where $ is a diffeomorphism near I x B 
such that 

$(i,0) = (t,0), 
®(0,y) = F Q - 1 (q(y)), 
$(t, y) = (t, y) for t > c with suitable c > 0. 

Write the components of q = Fq 1 o q as Taylor series 

qi(y) = q j (0) + V,f (0) • y + W{y)y ■ y 

where Hi are smooth matrices, and j = 0, . . . , n — 1 (t is the 0th 
variable). The properties of q imply that 

<f(0) = 0, fyg°(0) = 0, dpq a (0)=5%. 

We look for $ in the form 

&(t,y) = f(t) + a j (t) ■ y + Ri(t,y)y ■ y 

for some smooth functions j\ vectors a J and matrices TV . The condi- 
tions for $ motivate the following choices: 

f°(t)=t, T(t) = 0, a°M) = 0, a a Jt)=5%. 
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We choose R?(t,y) to be a smooth matrix with R?(0,y) = Hi (if) and 
BP(t,y) = for t > c. Then £)$(£, 0) = Id, and Lemma 17.31 ensures 
that $ is a diffeomorphism near Iq x B after possible decreasing B. □ 

The next result gives the construction of (non-reflected) Gaussian 
beam quasimodes associated with a finite length geodesic segment that 
enters and exits the domain nontangentially. To prepare for the re- 
flected case, we also consider the possibility of prescribing the bound- 
ary values of the quasimode at least up to high order at a point. Recall 
that if / is a smooth function having a critical point at p, the Hessian 
of / at p is the quadratic form 

Hess p (/)(7)(0), 77(0)) = (/ o V )"(0) 

where rj is any smooth curve with rj(0) = p. 

Proposition 7.5. Let 7 : [0, L] — > M be any unit speed geodesic in 
(M,g) such that j(0),j(L) £ dM, 7(0) and j(L) are nontangential, 
and 7(i) £ M mt for < t < L. Let also X be a fixed real number. For 
any K > there is a family {v s ; s = r + i\,r > 1} in C°°(M) such 
that as t — >■ 00 

||(-A 9 - s 2 )v s \\ L 2 (M) = 0(t~ k ), \\v s \\ L 2 iM) = 0(1), 

H^slU 2 (aAf) = 0(1), 

and for any if £ C(M) 

(7.1) I \v T+iX \ 2 ip dV g ->■ / e~ 2Xt ip(-f(t))dt as r ->■ 00. 

Jm Jo 

Given any neighborhood 0/7QO, L]) one can arrange that each v s is 
supported in this neighborhood, and away from the points where 7 self- 
intersects one has v s = e ts& a where 6 and a are smooth complex func- 
tions with 

dQ(r/(t)) = i(t)\ a( 7 (t)) ^ for r large. 

7/7 does not self-intersect at 7(0), the Kth order jets ofQ\g M and a\gM 
can be prescribed freely at 7(0) except for the following restrictions: 
dQ(j(0)) = 7(0) b , the Hessian o/Thi(6|,9m) at 7(0) is positive definite, 
and 0(7(0)) ^ 0. 

Proof. We embed (M, g) in a compact manifold (M, g) without bound- 
ary and extend 7 as a unit speed geodesic in M. Choose e > so that 
j(t) lies in M \ M and does not self-intersect for t £ [— 2s, 0) U (L, L + 
2e\ . We will construct a Gaussian beam quasimode in a neighborhood 
of 7 ([-£,L + e]). 
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Fix a point po = j(to) on j([—e,L + e]) and let (t,y) be any local 
coordinates near p , defined in U = {(t,y) ; t G /, \y\ < 5} for some 
open interval / containing to, such that p$ corresponds to (to,0) and 
the geodesic near p is given by T = {(t, 0) ; t G /}. Write x = (t,y) 
where x± — t and (x 2 , . . . ,x m ) = y. We seek to find a quasimode v s 
concentrated near T, having the form 

v s = e lse a 

where s = r + iX, and and a are smooth complex functions near T 
with a supported in {\y\ < 5/2}. 
We compute 

(-A - s 2 )v s = f 

where 

/ = e ise {s 2 [{(dQ, dG) - l)a] - is[2(dQ, da) + (A0)a] - Aa). 
We first choose so that 

(7.2) (dG, dG) = 1 to Nth order on T. 

In fact we look for of the form = J2f=o ®j where 

e i (*,y) = E § #i/ 7 - 

|7|=J 7 ' 

We also write = fl , / /c + 5jv+i where 

I/9N 

Set ^' fc = for Z > iV + 2. 

With the understanding that j, k run from 1 to m and a, (5 run from 
2 to m, the main part of the argument will consist of finding suitable 
©o, 0i and 2 in the following form: 

O (t) real valued, 
©i(t) = ia{t)y a with i a (t) real valued, 

2 (t) = l -H aP {t)y^ 

where H(t) = (H a p(t)) is a complex symmetric matrix, H a/3 = Hp a , 
such that lm(H(t)) is positive definite for all t. We also write 

6(0 = d t G (t). 
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Since <9 t 6 = £1 and <9 Q 0i = £ a , we compute 

g^djQdkQ - 1 = g l \d t Q + + . . .)(<9 t o + dS, + . . .) 

+ 2g la (d t & + $61 + . . .)(da&l + d a & 2 + ...) 

+ g a ^{d a Qi + d a Q 2 + . . .Xfye! +d p G 2 + ...)-l 

= g j %ik + 2g ll Ud t e 1 + ...)+ 2g la z 1 (d a e 2 + ...)+ 2 g la Ud t ®i + ■■■) 

+ 2g a ^ a (d^2 + ■■■) + <7 n (d t 0i + • • O(^0i + • • •) 
+ 2g la (d t Q 1 + . . 0(9*02 + • • + 9 aP (d a Q 2 + . . 009/302 + • • " 1 
= 9 j %tk+2g 1 %(dt6i+. • .)+2g a %{d a e 2 +. . .)+g 11 (d t G 1 +. . .)(d&i+- . 
+ 2<? lc *(c> t 1 + . . .){d a B 2 + ...)+ g al3 (d a Q 2 + . . 0(^02 + . . - 1- 

Writing g^ k = g 3 h + g{ k + . . . and grouping like powers of y, we obtain 

(7.3) gi%ed k e - 1 = [gi%^ k - 1] 

+ [g{%ik + 2g 1 %i^ + 2gS%H a pyP} 
+ {92 + ■ ■ -Mk + 2g 1 %(d t Q 2 + ...) + 2{g\ k + . . 0^(^01 + . . 
+2g%%(d a e 3 +. . .)+2(gf+. . O&(d Q 2 +. . O+(? 11 («9 t 1 +. . O(d t 0i+. . 

+ 2g la (d t Q 1 + . . .){d a e 2 + ...)+ g^(d a Q 2 + . . O(^0 2 + . . 0- 

We can make the two expressions in brackets to vanish by choosing 
£(t) to be part of the solution (x(t),£(t)) of the cogeodesic flow with 
Hamiltonian h(x,£) = \g 3k {x)^ k , 

x i (t) = dt.h(x(tU(t)), 
Ut) = -d XJ h(x(t),t(t)). 

There is a unique solution with x(t ) = p and £(to) = 7(to) b (here we 
raise and lower indices with respect to the metric g). It follows that 
x(t) is the unit speed geodesic t h-> (t, 0), and ^(t) = x 3 \t). Then 
g^ijik = 1, and with our choice of coordinates = 1 and £ Q = so 
that also 

g]% = l, gf£ k = 0. 

We further have 

ipv" = -\d x ,g jk {tM^ = -\g{%ik- 

Noting that d\ has unit length, we have 

Z 1 = g lk (t,0)Z k = l. 
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Since £ Q = g a k(t, 0)£ fc = g a i(t, 0), we can therefore choose 

e (t) = t, 

Q 1 (t)=g al (t,0)y a . 

Using these choices and the facts above, in (17. 3ft the expressions in 
brackets will indeed vanish and one obtains 

s^-e&e-i = (gi k +. . .M k +2(d t e 2 +. . .)+2( 5l lfe +. . .)&(d t e 1 +. . .) 
+ 2( g f + . . .)&(a«e 2 + ...) + <7 n (atei + . . + . . .) 
+ 2 9 la (d t e 1 + . . .)(d a e 2 + ...) + g a ?{d a e 2 + . . .)(fye 2 + . . .) 

= d%Zk + 2d t e 2 + 2 9 f£ k d a e 2 + 2g 1 a d t e 1 d a e 2 + g fd a e 2 d^e 2 

+ 2gl%d t Q 1 + g 1 1 (d t Q 1 ) 



+ Yl + 2d t e p + 2g?£ k d a e p + 2 9o - a d t e 1 d a e p + 2 g fd a e 2 d^e p 

p=3 

p—1 p—1 P — 2 

+ 2Y,9l- l ikd t Q l + 2Y^9t: i+l ikd a Q l + Y^9] 1 E d ^rd t Q s 

1=1 1 = 2 1=0 r+s=p-l 

l<r,s<p 

p-2 p-2 

+I> lQ E ^ e Ae,+E^ E d a e r d^e s ]+o(\y\ N+1 ). 

1=0 r+s=p-l+l 1=0 r+s=p-l+2 

l<r<p 2<r,s<p 
2<S<p 

We want to choose G 2 so that the first term in brackets vanishes. 
Recalling that we are looking for 6 2 in the form Q 2 (t, y) = ^H a/3 {t)y a y^ 
where H(t) is a smooth complex symmetric matrix, it follows that H 
should satisfy the matrix equation 

iZ^l/Y + 2gJ%H^y p + 2g^d t Q l H 1 ^ + gf H ia H 5p y a y p = F aP y a y? 

where F(t) is a real valued smooth symmetric matrix. This can be 
further written as the matrix Riccati equation 

H + BH + HB t + HCH = F 

where B(t) and C(t) are real smooth matrices and C is symmetric. 
More precisely, since g{ = d a g^ h (t, 0)y a we have 

(7.4) Bl = d a g^ k {tMk + gVL Cn 5 = gl 5 . 

Choosing H(to) = Ho where Hq is a complex symmetric matrix with 
lm(H ) positive definite, it follows that the Riccati equation has a 
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unique smooth complex symmetric solution H(t) with lm(H(t)) pos- 
itive definite |KKL01j . This completes the construction of 02- From 
the lower order terms we can find G3, . . . , 0jv successively by solving 
linear first order ODEs on T with prescribed initial conditions at to. In 
this way we obtain a smooth satisfying ( 17. 2p . 
The next step is to find a such that 

s[2(de, da) + (A8)o] - iAa = to iVth order on T. 

We look for a in the form 

a = r * {a + s a-i + ... + s a_ N )x{y/oj 

where x is a smooth function with ^ = 1 for |y| < 1/4 and \ = for 
\y\ > 1/2. Writing r\ = A0, it is sufficient to determine aj so that 

2(dQ, da ) + r]a = to Nth order on T, 

2(dQ, da_i) + rja^i — iAao = to Nth order on T, 

2(dQ, da_x) + V a -N — iAa-(N-i) =0 to iVth order on T. 

Consider ao = aoo + • • • + cion where aoj(t, y) is a polynomial of order 
j in y, and similarly let r\ = r] + . . . + r/ N . We compute 

2(dQ, da ) + 7]a = 2{g 1 1 + . . .)(<9 t Q + . . .){d t a 00 + . . .) 

+ 2(g 1 a + ...)(d t Q + ...)(d a a 01 + ...) 
+ 2(g 1 a + ...)(d a e 1 + ...)(d t a 00 + ...) 
+ 2(gf + ...)(dpG 1 + ...)(d a a 01 + ...) + (r ]o + r] 1 + ...){a QO + a i + ...). 

Recalling that d t Qo = £1 = 1, 9 a 9i = £ Q where <7 3 £j = 1 and <?o 3 £j = 0> 
we obtain 

2(d9, da ) + V a = 2 [gl% + ^ Q u (9 t 0i + ...) + (ft 11 + . . .){d t &o + . . .) 

+ gl^a + gl a {d a ®2 + •••) + (gl a + ■ • -X^i + •••)] (9 t a 00 + ...) 

+ 2[gl a (d t e 1 + ...) + (g{ a + . . .)(^0 O + ...)+ gf(dp<d 2 + . . .) 
+ + . . .)(fy©i + ...)] {d a a 01 + ...) + (ijo + Vi + ■ • OKo + a i + . . .) 

iV 

= [2d t a 00 + r] a 00 ] + ^2 hd t a 0p +qfy^d a a 0p +r]oa 0p +F p + 0(\y\ N+1 ) 
P =i 

where q p {t) are smooth functions only depending on g and 0, and 
F p (t,y) is a polynomial of degree p in y that only depends on g, 0, 77 
and aoo, • • • > a o,p-i- 
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We want to choose aoo so that the first term in brackets vanishes, 
that is, 

^a o + T^oaoo = 0. 

This has the solution 

a 00 (t) = c e-^<> (s)ds , a 00 (t ) = c . 

We obtain a i, . . . , a Ar successively by solving linear first order ODEs 
with prescribed initial conditions at to- The functions ai, . . . , ajv may 
be determined in a similar way so that the required equations are sat- 
isfied to iVth order on T. This completes the construction of a. 

To review what has been achieved so far, we have constructed a 
function v s = e tse a in U where 

0(t, y) = t + Ut)y a + l -H a p{t)y a y? + 9, 

a(t,y) = r 2r V L (a + s _1 a_i + . . . + s~ N a_ N )x(y/S'), 

a o (t,0) = c o e-^> ois)ds . 

Here = 0(|y| 3 ) and and each aj are independent of r. Also, 
/ = (—A — s 2 )v s has the form 

/ = e isQ r^ (s 2 h 2 a + SIH + ... + s-^-^/i-^-i) +is~ N A(a- NX (y/6'))) 

where for each j one has hj = to iVth order on T. We also note that 
de( 7 (t)) = i(tf and Hess 7(to) (lm(0| {t=to} )) = lm(H (*„)). 
To prove the norm estimates for v s in U, note that 

\ e is °\ = e -ARce e -rIme = e -\t e -\rl m {H{t))y-y e -\0{\y\) e -rO(\y\ A ) _ 

Here lm(H(t))y ■ y > c\y\ 2 for (t,y) £ U where c > depends on H Q , 
g and I. By decreasing 5' if necessary, this shows the following bound 
when t in a fixed compact set: 

\v s (t,y)\<r^e-^ 2 x(y/5'). 

Integrating the square of this over U we get, as r — > oo, 

INU 2 (t/) ^ llr^e - ^! 2 !!^^ = 0(1). 

Similarly we have 

||(-A - <?)v.\\v m < ||r a ¥ 1 e-i"l''l 5 (r 2 | ! /r +1 + T-")\\ mu) 

= 0(r^). 

The norm estimates for v s in U follow upon replacing iV by IK + 3. 

For the L 2 (dM) estimate, if U contains a boundary point x = 
(to,0) G dM, then by assumption ^| xo is transversal to dM. If p 
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is a boundary defining function for M, so dM is given as the zero set 
p(t,y) = near xq and Vp = — v on dM, then |f(^o) 7^ and by the 
implicit function theorem there is a smooth function y \-± t(y) near 
such that dM is given by {(t(y), y) ; \y\ < r } near x . The bound for 
v s given above implies that for 8' small 

\\ v s\\h {dM nu) = / l^(*(2/)>2/)l^(l/) 

]j/|<»"0 

< f T^e^\y\ 2 d y = 0{l). 

J\y\<r 

At this point we can construct the quasimode v s in M from the 
corresponding quasimodes defined on small pieces. Let "f([—e,L + e)) 
be covered by open sets U^°\ . . . , [/( iV+1 ) as in Lemma I7.4[ and note 
that each corresponds to Ij x B(0,5) in the (t, y) coordinates. 
Suppose first that 7 does not self-intersect at time t = 0. We find a 
quasimode 

t,(o) = e «e {0) a (°> in f/ ( ' ^ by the above procedure, with some 
fixed initial conditions at t — for the ODEs determining 0(°) and a^°\ 
Choose some t' with 7(^0) G t/( ) fl tA 1 ), and construct a quasimode 
_ e ise 1 ) a (i) j n choosing the initial conditions for the ODEs 

for and at t' Q to be the corresponding values of G*- -* and at 
t' . Continuing in this way we obtain v^ 2 \ . . . , v ( N+1 \ If 7 self-intersects 
at t — 0, we start the construction from fixing initial conditions for 
the ODEs at t = 0, and find v (0) by going backward. 

Let {Xj(t)} be a partition of unity near [— e, L + e] corresponding to 
the cover {Ij}, let Xj{t,y) = Xj{t) 011 and define 

JV+l 
j=0 

Then v s is smooth in M and it is supported in a small neighborhood 
of j([—e,L + e]). The important point is that since the ODEs for 
the phase functions and amplitudes have the same initial data in U^> 
and in U^ +l \ and since the local coordinates (fj and v^j+i coincide on 
<fj 1 ((l j nlj +1 ) xB), one actually has v {j) = v {j+1) in (pj 1 ((l j nl j+1 ) xfl). 
Letting pi, . . . ,p R be the points where 7 intersects itself, we choose an 
open cover of supp(t> s ) fl M, 

supp(^) n M c (u* =1 K) u (uf+o^H^-o u w jA )) , 
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where V r are small neighborhoods the points p r and Wj fi , Wj ; x C 
such that 

v s\v r = vi3) > 

l(tj)=Pr 

Vs\w jtl =V^ +l \ 

Since v s is a finite sum of the v^> in each case, the L 2 (M) bounds 
for v s and (—A — s 2 )v s and the L 2 (dM) bounds for v s follow from 
corresponding bounds for the The form of v s near points where 
7 does not self-intersect and the possibility to prescribe the Kt\i order 
jets of Q\dM and cl\qm at 7(0) follow from the construction and Lemma 

m 

To conclude the proof, using a partition of unity it is enough to verify 
the limit (17. ip for any ip supported in one of the sets V r P\M or H^/flM. 
Further, we can choose the sets V r to be so small that the real part of 
dQti) — is nonvanishing near V r if 7(i?) = 7(tfc) = Pr but j 7^ k. 

This follows since 

Re ( rf0 O) _ d<d^)( Pr ) = 7 (t i ) b - 7 (t fc ) b ^ 0. 

Here we may need to decrease 5 so that we still have an open cover. 

Consider first the case where if) G C c (Wj t i fl M). Here the support 
of ip may reach <9M, and we extend ip by zero outside of Wjj fl M. 
Suppose that v s = e tsB a where 6 = t + ^ a y a + ^H(t)y ■ y + 0(\y\ 3 ) and 

a = r Z2 r i (a + 0(T- l ))x(y / '5'), and let p = \g\ l l 2 . Then 
\v T+i x\ 2 ip dV g 



M 

00 



r 2At e -rM^))^ e TO(^ 



OO 



e -2At / e -lm(H(t))yy e T-^0(\yf) e r-^0(\y\) x 



(\a (t, T" l/2 y)\ 2 + 0{T~ l )) X (y/T l l 2 5') 2 ^ r'^pit, r~ l ' 2 y) dt dy. 

Since lm(H(t)) is positive definite and 5' is sufficiently small, the term 
e -im(H(t))y-y dominates the other exponentials and one obtains 



>3C 



lim / \v T+iX \ 2 ip dV, 



M 

2Xt 



e ~~ { I e -iMH{t))vv dy ) \a (t,0)\ 2 ip{t,0)p{t,0)dt. 
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Evaluating the integral over y gives 

lira / |^|y dV 9 = C m f e-^ '"f ' °™\%(t, 0) dt. 

We will prove below that 

{7 |a (t,0)|V(t, 0) 

(7.5) — ; = const. 

y/det lm(H(t)) 

The limit ( 17. ip will follow upon dividing the family {v s } by a suitable 
constant. 

If ip £ C c (V r Pi M) (again supp(-0) may extend up to dM), we have 

v s\v r = V ^> 

7(*j)=Pr 

so that on V r 

\ v f= \v U) \ 2 + J2 v(j) ^ m - 

y(tj)=p T l{tj)=l{tk)=Prd¥=k 

We arranged earlier that Re(c?0^ — dQ^) is nonvanishing near V r if 
7(^') = 7(*fc) — Vr but j ^ k. Thus the cross terms give rise to terms 
of the form 



v (J) v (k)ip dV = / e iTip w u) w^ipdV 
VrtlM Jv r nM 

where = Re(0^ — 0( fc )) has nonvanishing gradient in V r , and = 

e i,Im(0( i )) e -ARe(*(O) a (O > Wg wigh ^ prQve th&t 

(7.6) lim ! e* v w®wWif>dV = 0, 3^ k, 

T ^°° Jv r nM 

showing that the cross terms vanish in the limit and the previous com- 
putation for |t>^| 2 shows the limit (17. ip also when ip is supported in 
some V r fl M. To show (17.61) . let e > 0, and decompose ip — 4>i + ^2 
where ip x G C C °°(K- n M) and ||^||i»(v P nM) < £■ Then 



e ^ w 0') w (fc)^ 2 < \\ w {j) \\ L 2 \\w {k) \\ L 2 1| 

v r nM 

since ||u;0||xa < H^Hl 2 i$ 1- F° r the smooth part ipi, we employ a 
non-stationary phase argument and integrate by parts using that 

jrv = —L(e lTV ), Lw = (\d(p\- 2 d(p, dw). 

IT 
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This gives 

e^w^wMfadV = I ^^ (]) #)^ dS 
v r nM JdM tT\dtp\ 2 

Since \\v^\\ l 2 {8m) — 0(1), the boundary term can be made arbitrarily 
small as r — > oo. As for the last term, the worst behavior is when 
the transpose L l acts on e lsIm ( e(!) ) ? anc l these terms have bounds of the 
form 

|||d(im(e^))|^|| i2 ||^ fe )|U 2 |^i|U- 

Here |d(lm(0^))| < \y\ if (t, y) are coordinates along the geodesic 
segment corresponding to and the computation above for H^^Hl 2 
shows that 

IlldOmCe^i^ii^ii^ii^ii^lli- < r- 1 ' 2 . 

This finishes the proof of (17.61) and also of (17. ip . 
It remains to show (17. 5p . We have 

|a (t,0)| 2 pM) = \co\ 2 e-^ M(s) %(t,0)\ 1/2 - 
Note that r] (t) is given by 
%(*) = A9(t,0) 

= g u d 2 Q + 2g la d ta Q 1 + g^d aP e 2 + d jg jl d t Q + d^d^ 
+ ±d j (\og\g\)(gi 1 d t e + g^d a e 1 ) 

= 2g la i a + g^Hap + ( V% + ^(logM)^ 

The conditions gi% = 5{ and g la £ a = g lk £ k = -(d t g lk )£ k at (i,0), 
together with the general fact that d t ()og\g\) = —gjkdtg^, imply that 

Vo(t) = g la L + g aP H afi + (d a g ak )i k + l -d t {\og\g\). 

Recalling the definition of the B and C matrices in (I7.4p . this says 
precisely that 

r ]0 (t) = B2 + C^H ya + ^d t (\og\g\) 

= tr(B(t) + C(t)H(t)) + ±d t Qog\g\). 
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Consequently | o (t, 0)\ 2 p (t, 0) = c' Q e~ St ^{ s )+c{s)^ms)))ds Qn the 
other hand, by |KKL01[ Lemma 2.58] solutions of the matrix Riccati 
equation have the property that 

detlm(tf(t)) = detIm(i/(to))e- 2 ^ to ( B W +c ^ Re ^^)) ds . 

This proves the result. □ 

The proof of Proposition [771] now follows rather quickly from the way 
we have set up the previous result. 



Proof of Proposition 7.1. Let 7 : [0, L] — > M be a nontangential broken 



ray with endpoints on E, and let < t\ < . . . < t N < L be the times of 
reflection. Let vi ^ be a Gaussian beam quasimode as in Proposition 17. 51 
associated with the geodesic 7|[o,ti] • We will construct another Gaussian 

beam quasimode vi associated with 7|[t ll t 2 ] such that v& — Vs \qm w iU 
be small near 7(^1 ). 

In fact, by Proposition 17.51 we have vi = e lse{j) near 7(^1), and 
we can choose the Kth order jet of Q^'\dM at 7(^1) to be equal to that 
of 0(°)| qm with the following exception: we always have 



rf(© (0) )U)=7(ti-) 



rf(@ (1) )U)=7(*i+) b , 



It follows that 

d(e {0) \ d M)U(t l) = i(t 1 -) 
d(e% M )| 7(M = 7 (*i+) 



where we have taken the projections to the cotangent space of dM 
at 7(£i). But by the rule that the angle of incidence equals angle of 
reflection, it holds that 7(^1— ) tan = j(ti+) ta n- Thus actually the Kth. 
order jets of Q^°'\dM and Q^\dM coincide at 7(^1), and by Proposition 
17.51 we can also arrange that the i^th order jets of aS Q ^\dM and cl^\qm 
coincide at j(ti). 

Write f s = — Vs\qm, and let (t, y) be coordinates near 7(^1) 
such that dM is parametrized by y H- (ti,y) and 7(^1) corresponds to 
(ti,0). Recall that Vs are supported in small tubular neighborhoods 
of the corresponding geodesic segments. By the above considerations 
and the construction of 0^' and and dropping the variable t% from 
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the notations, the restrictions of 0^ and to dM satisfy 

eW(y) = 9(y) + sW(y), a^(y) = a(y) + 6^(y) 

where is a polynomial of order K, a = r I2 4 1 ax(y/5') where a is 
a polynomial of order K, and |S^(y)| < C\y\ K+1 and |6^(y)| < 
Cr I2 ^|y| is ' +1 x(y/5') on supp(x( • /$'), where \ is a cutoff function and 
5' is a constant independent of r that can be chosen as small as we 
want (these initially depend on j, but since there are finitely many 
reflections we can choose them independent of j). Here and are 
independent of r, and a and 6^0 are mildly r-dependent and satisfy 
uniform bounds with respect to r. Then 

f s = e ise ((e isE(0) - e lsSW )a + e^V* - e^'V*). 

We have 

and consequently near y = 

|e isS(0) -e isS(1) |<Cr|y|^ +1 e c ^l K+1 . 

Thus, near y = 

|/ s (y)| < Cr^e-Me) r|yr+ i e c.|^ +lx(y/50 _ 

Using that the Hessian of lm(0) at is positive definite and choosing 
5' sufficiently small, we have 

\f s {y)\<Cr^e-^y\\\y\ K+1 X {y/5'). 

Integrating the square of |/ s | over W n ~ l and changing y to r _1 / 2 y, we 
obtain 

\\fs\\L HRl) = 0(r' E ^) 

where Ri is a small neighborhood of 7(^1) on dM containing the set of 
interest. 

Repeating this construction for the other points of reflection, we end 
up with a quasimode 

N 

v s = £(-l)'"t;?> 

3=0 

that is supported in a small neighborhood of the broken ray 7. Since 
all points of reflection are distinct, we can arrange that the quasimode 
satisfies 
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It also satisfies 

||(-A - s 2 )v s \\ L 2 (M) = 0(t- k ), \\v s \\ L 2 {m) = 0(1). 

Replacing K by 2K + 1, we have proved all the other statements in 
the proposition except for the expression of the limit measure. To 
do this, we consider the finitely many points where the full broken 
ray 7 self-intersects or reflects, and decompose the terms as in 
the proof of Proposition 17.51 to parts living in small neighborhoods of 
the self-intersection and reflection points and parts away from these 
points. Now all self-intersection points are in the interior or on E and 
all self-intersections must be transversal, and also all reflections are 
transversal. Consequently, when forming \v s \ 2 , the cross terms arising 
from different parts living near the same self-intersection or reflection 
point contribute an o(l) term by non- stationary phase as in the proof 
of Proposition 17.51 Thus the limit measure of |f s | 2 (iV^ is indeed the 
measure e _2A< 5 7 , where <5 7 is the delta function of the broken ray 7. □ 

8. Recovering the broken ray transform 

In this section we give the proof of Theorem 12.41 concerning the re- 
covery of integrals over broken rays. 

Proof of Theorem \2.4\ The proof is very similar to the proof of Theo- 



rem 12.11 except that we use reflected Gaussian beam quasimodes in- 
stead of WKB type quasimodes. Let 7 : [0, L] — > M be a nontangential 
broken ray with endpoints on E, and let A > 0. Let also g = e © go 
and ijj = c(qj — c 1 ^ A 9 (c _I1 ^)). Consider the complex frequency 

s = r + i\ 

where r > will be large. We look for solutions 

u x = e- SXl (v s (x') +n), 
u 2 = e SXl (v s (x')+r 2 ), 

of the equations (— + qi)u\ = 0, (— A^ + q2)u2 = in M. Here 
v s G C°°(M ) is the quasimode constructed in Proposition 17.11 that 
concentrates near the given broken ray 7 and is small on OMq \ E. 
Since Ag = df + A 90 , the function u\ is a solution if and only if 

e SXl {-A g + gi)(e- Wil ri) = -(-A go + ~q x - s 2 )v s (x') in M. 

We look for a solution in the form 77 = e tXxi r[ where r[ satisfies 

e TX1 (-Ag + q 1 )(e- TX W' 1 ) = f in M 

with 

/ = -e-^i-A^ + q,- s 2 )v s (x'). 



f- 



50 CARLOS KENIG AND MIKKO SALO 

To arrange that Ui\r 4 = 0, fix some small 5 > 0, let S± and Sq be the 
sets in Proposition 14. 31 with Carleman weight <p(x) = —X\, and consider 
the boundary condition 

e mIs_us — e J- 

where 

-e~ lXxi v s (x'), on Tj, 

0, on (5_U5o)\r,. 

For any fixed K > 0, by Proposition 17.11 and by the condition that 
r.j C R x (dMo\E) we may assume that the following bounds are valid: 

||/|| L2(M ) = 0(1), ||/_|U 2( 5_) = 0, ||/-|U 2( 5 ) = 0(r- K ). 

It follows from Proposition 14. 31 that there is a solution r[ satisfying the 
above boundary condition and having the estimate 

IKHl2(a/) = 0(t~ 1 ). 

Choosing r[ as described above and choosing r\ = e lAxi r^, we have 
produced a solution U\ G H&.(M) of the equation (— A^ + qi)u\ = in 
M, having the form 

fii = e- SXl {v s (x') +n) 

and satisfying 

supp(tti| aM ) C dM + U <9M_ U T a 

and 1 1 T"a 1 1 z,a (at) = 0(r _1 ) as r — > oo. Repeating this construction for 
the Carleman weight </ , (^) = a^i, we obtain a solution u 2 G Ha-(M) of 
the equation (— + ^2)^2 = in M, having the form 

U2 = e SXl (v s (x') + r 2 ) 
satisfying the same support condition and bound for ||r"2 ||x 2 CJWf) ■ 

n _2 

Writing Uj = c~^Uj, Lemma 16.21 shows that Uj G H& g (M) are 
solutions of (— A g + q\)u\ = and (— A g + q2)u 2 = in M. Then 
Proposition 16.11 implies that 

(qi - q 2 )uiU^dV g = 0. 

M 

We extend qi — q 2 by zero to R x Mq. Inserting the expressions for Uj, 
and using that dV g = c n ^ 2 dx\ dV go (x'), we obtain 



Mo 



(qi-q2)ce- 2lXxi ( y \v s (x')\ 2 + v s r 2 - + — s r 1 +r 1 rl) dx x dV go (x') = 0. 



Since ||rj||L3(M) = 0(r x ) as r — >■ 00, Proposition 17.11 implies that 



L 



-2\t, 



c{q 1 -q 2 )y{2\ 1 {t))dt = 0. 







CALDERON PROBLEM WITH PARTIAL DATA 



51 



This concludes the proof. 



□ 
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